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Abstract 


An  analytical  study,  using  the  STAGSC-1  computer  code,  was 
conducted  on  a  Kevlar/Polyester  composite  shell  of  general  shape  with 
internal  asymmetrical  pressure  loading.  Experimentation  was 
conducted  on  the  shell  -by  feho-Aii^  force  Wright  Aeronautics 
Labroratery ,  Wright -Pafterjen  Afbi — Ghi^to  find  displacements  in  the 
shell  due  to  this  internal  pressure  loading. 

Modeling  of  the  composite  structure  was  done  in  increments 
whereby  each  change  in  the  finite  element  model  better  approximated 
the  actual  shell.  Nonlinear  computer  runs  were  done  at  each  model 
increment  to  compare  against  experimental  results.  Linear  computer 
runs  were  also  completed  for  comparison  purposes. 

It  was  found  that  accurate  modeling  of  the  shell  to  include 
thickness  variations,  boundary  conditions,  and  materials  is  essential 
to  obtain  reasonable  results.  Also,  the  incorporation  of  the 
nonlinear  analysis  leads  to  displacement  results  that  are  within  15% 
of  experimentation.  Linear  results  from  the  same  model  are  in  error 
by  over  75%  due  to  large  displacements  in  the  shell.  ! 

r  'vVv^  c  K  ^  ^  H  t 

{v  a  '"V  iVi-U  .  'I  . 
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NONLINEAR  FINITE  ELEMENT  ANALYSIS  OF  A  GENERAL 
COMPOSITE  SHELL 

1 .  Introduction 

In  this  thesis,  a  large  laminated  composite  shell  is  analyzed 
using  a  finite  element  program.  Since,  under  loading,  large 
displacements  occur  in  the  shell,  a  nonlinear  approach  is  used  in  the 
analysis.  Shell  research,  in  general,  is  important  since  this  basic 
component  is  in  many  aerospace  vehicles.  For  example,  the  fuselage, 
wings,  and  canopy  on  an  aircraft  or  the  skin  on  a  missile. 
Furthermore,  there  is  an  increased  use  of  composite  materials  in  many 
of  these  shell  components  as  in  the  use  of  composite  wings  on 
aircraft  or  the  casing  of  a  solid  propellant  rocket  motor.  This  is 
primarily  due  to  the  ability  of  a  designer  to  tailor  the  material 
properties  in  a  structure  to  withstand  the  load  environment.  This 
material  tailoring  leads  to,  in  general,  a  lighter  structure  than 
similar  metal  structure  thereby  reducing  the  overall  weight  of  a 
vehicle;  a  critical  factor  in  aerospace  design. 

The  basic  area  of  work  on  shells  is  very  broad.  A  good  portion 
of  this  work  deals  with  geometrically  linear  problem  solving.  Ugural 
[1]  shows  a  basic  layout  of  linear  shell  theory.  Other  works  such  as 
Saada  [2]  and  Sanders  [3]  deal  with  shells  from  a  nonlinear  point  of 
view. 

Fewer  references  consider  composite  shells  modeled  nonlinearly 
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and  most  of  tnese  deal  with  a  standard  geometric  shape  such  as  a 
cylinder.  Bauld  [4]  and  Dennis  [5]  work  with  shells  of  this  type. 
Dennis  [5]  also  has  a  very  good  introduction  that  outlines,  in 
detail,  other  works  (over  100)  in  shell  analysis. 

With  the  ever  increasing  capability  of  the  modern  computer,  many 
shells  are  analyzed  numerically.  Most  of  the  recent  computer  work  is 
finite  element  based.  Cook  [6]  has  a  very  well  written  section  on 
nonlinear  solution  techniques  while  the  Structural  Analysis  for 
General  Shells  (STAGS)  theory  manual  [7]  and  Zienkiewicz  [8]  have 
sections  that  deal  specifically  with  shells  and  geometrically 
nonlinear  finite  element  solutions. 

Other  references  are  concerned  with  the  experimentation  on 
composite  shell  structures.  This  is  of  particular  importance  in 
order  to  validate  closed  form  and  finite  element  solutions. 

References  in  this  area  include  works  by  Knight  and  Starnes  [9] , 
Tisler  [10],  and  Lee  [11]. 

This  thesis  addresses  all  of  these  aspects  of  composite  shells 
and  includes  another  aspect  that  is  not  well  documented.  This  other 
aspect  is  the  comparison  of  a  large  scale  nonlinear  analysis  with  the 
experimental  testing  of  a  general  composite  shell.  There  is  very 
little  documentation  on  shells  of  the  size  (approximately  18  ft. 
long,  3  ft.  wide,  and  3  ft.  high)  and  shape  (not  a  standard  geometric 
shape)  considered  in  this  thesis. 
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1 . 1  Objectives 

The  major  purpose  of  this  thesis  was  to  determine  If  a  nonlinear 
finite  element  analysis  of  a  shell  produces  better  results  than  a 
linear  analysis  of  the  same  shell.  The  program  used  to  do  this 
analysis  is  the  Structural  Analysis  for  General  Shells  (STAGSC-1) 
computer  program.  The  shell  is  loaded  with  an  internal  asymmetric 
pressure  load.  The  second  purpose  was  to  investigate,  the  effect  of 
modeling  the  shell  as  accurately  as  possible  to  see  how  changing 
different  model  parameters  effect  the  overall  shell  response.  A  third 
purpose  was  to  compare  finite  element  results  to  experimental  values. 
The  final  purpose  was  to  carry  out  a  stress  analysis  to  check  for 
composite  failure.  The  results  will  be  used  to  better  understand  the 
effects  of  large  displacements  on  a  finite  element  analysis. 


1 . 2  Scope 

A  total  of  eight  models  were  run  on  the  shell  to  find 
displacements  due  to  an  internal  asymmetric  pressure  load.  The 
models  were  developed  incrementally  by  starting  with  a  constant 
thickness  isotropic  shell  with  clamped  boundary  conditions.  Each 
model  increment  changed  the  basic  model  to  better  approximate  the 
actual  composite  shell.  These  changes  included  modeling  the 
thickness  variation,  modeling  the  composite  material,  and  modeling 
the  boundary  conditions  which  are  not  a  true  clamped  condition. 
During  these  increments  the  model  was  run  nonl inear ly  using  STAGSC-1 
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with  linear  runs  for  comparison  purposes. 

Next,  the  results  were  compared  to  experimental  results.  The 
experimental  results  were  obtained  by  the  Wright  Aeronautics 
Laboratory,  Wright -Patterson  AFB,  Ohio. 

Finally,  a  Tsai-Wu  failure  analysis  was  carried  out  on  the  shell 
to  determine  if  ply  failure  occurred  in  the  composite. 
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2.  THEORY 


The  Structural  Analysis  for  General  Shells  (STAGS)  computer 
program,  developed  by  Lockheed  Palo  Alto  Research  Lab,  was  first 
operational  in  1967.  Its  primary  purpose,  as  evident  in  its  name,  is 
the  analysis  of  thin  shelled  structures.  From  1967-1976  the  program 
was  based  on  the  finite  difference  method.  In  1979  a  new  version  of 
STAGS,  STAGSC-1,  was  released  and  is  based  entirely  on  the  finite 
element  method  [12].  The  1986  STAGSC-1  VAX  Computer  version  is  used 
in  this  thesis.  Of  the  many  capabilities  in  the  STAGS  program,  the 
primary  one  used  herein  is  the  nonlinear  static  analysis.  The  basic 
finite  element  used  in  STAGSC-1  is  the  flat  plate  element  that  facets 
the  shell  surface  to  approximate  its  curvature .  A  look  at  the  theory 
for  this  flat  plate  element  and  the  theory  used  in  STAGS  is 
presented  in  the  following  sections.  VJhile  this  theory  may  not  be 
the  exact  theory  used  in  STAGSC-1,  it  is  an  attempt  through  research 
of  documentation  to  give  the  reader  some  insight  to  STAGS'  internal 
make-up  to  better  understand  the  code. 

Also,  a  review  of  Sanders'  nonlinear  shell  theory  is  presented 
to  better  understand  the  kinematic  relations  needed  in  a  plate 
element  to  approximate  a  shell  surface.  This  does  not  imply  the  use 
of  Sanders'  nonlinear  kinematic  relations  in  STAGS;  its  is  just  a  way 
of  showing  important  relations  in  a  true  shell  formulation  that  need 
to  be  incorporated  in  a  plate  type  element. 
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2 . 1  Nonlinear  Shell  Theory 


A  review  of  Sanders’  nonlinear  shell  theory  is  presented  in 
order  to  understand  the  nonlinear  strain  displacement  relations  in  a 
shell.  This  presentation  will,  in  the  subsequent  section,  provide  an 
insight  to  the  theory  used  in  STAGSC-1. 

The  classical  thin  shell  theory  derived  by  Sanders  assumes  that 
the  shell  is  thin,  middle  surface  strains  and  rotations  are  small, 
and  displacements  away  from  the  midsurface  are  restricted  by  the 
Kirchhof f -Love  hypotheses  [3].  With  these  assumptions  in  mind  the 
equations  for  the  midsurface  and  bending  strains,  respectively,  are 

o  X.  x,y  W  1,2.  1,2 

£  .  -  +  - i.  v  +  -  +  t:  <f>  +  4> 

X  a  QiQ  R2x2 

X  X  y  X 


o  y  .  y.x  .w,l,2  1,2 

e  -  — -  +  — —  u  +  -  +  x  4  ■  4 

y  a  aa  R2y2 

y  X  y  y 


V,  u. 


V  +  4  4 


(2.1) 


where  u  and  v  denote  displacements  tangent  to  the  midsurface  and  w 

denotes  displacement  perpendicular  to  the  shells  midsurface,  R  and 

R  are  the  principal  radii  of  curvature,  4  >  4  >  and  4  are  in-plane 
y  X  y 
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rotational  terns,  and  a  and  a  are  metric  coefficients, 

*  y 

rotational  terms  are  given  by 

w, 


The 


X 


a  R 

X  X 


W, 

_ I  +  X 

a  R 
y  y 


(2.3) 


*  -  2^  [<V>  ■  •  •,] 

X  y 


The  metric  coefficients,  written  in  a  form  similar  to  Saada  [2],  are 


-  JE  (  1  .  ^  ) 

-  JC  (  1  .  ^  ] 


(2.4) 


where  is  Che  distance  from  the  midsurface  and  E  and  G  are 
functions  of  the  shell  geometry  (not  material  constants) . 

In  order  to  simplify  the  presentation  and  get  a  better  feel  for 
the  strain  equations,  a  thin  cylindrical  shell  will  be  used  (see 
Figure  2.1).  The  reference  surface  for  the  shell  is  the  midplane. 
For  this  cylindrical  shell  we  let  R  -««>,  R -R,  f  -z ,  and  assume  that 

X  y  3 

z<  R.  For  this  particular  case  Sanders'  equations,  (2.1)  and  (2.2), 
reduce  to  [4] 

£°  -  u,  +  1/2  +  1/2  / 

XX  X 


£°  -  V,  +  ^  +  1/2  4^  -  1/2  / 

y  y  R  y 


(2.5) 


and 


7  -  V,  +  u,  +  <f>  <f> 

xy  X  y  X  y 


K  —  0 
X  x,x 


(2.6) 
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Figure  2.1  Cylindrical  Shell  Section  Showing  Coordinate 
and  Displacement  Directions  [4] 


Now,  equations  (2.5)  and  (2.6)  can  be  written  with  the  Kirchhoff -Love 
hypothesis  to  obtain  the  full  strain  expression; 


From  equations  (2.5),  (2.6),  and  (2.7)  it  is  much  easier  to 


understand  the  strain  displacement  relations  for  a  simple  shell  while 
still  appreciating  the  complexity  of  the  formulation. 


2.2  STAGSC-1  Theory 

2.2.1  Strain  Displacement  Relations 

As  stated  earlier,  STAGS  uses  flat  plates  to  approximate  the 

curved  surface  of  a  shell.  These  plate  elements  are  thin,  therefore 

a  state  of  plane  stress  can  be  assumed  with  7  ,7  ,  £  ,  and  a 

xz  yz  z  z 

equal  to  zero,  and  the  in-plane  displacements,  u  and  v,  as  well  as 
the  normal  displacement,  w,  functionally  depending  on  only  two  space 
variables  [7].  As  in  Sanders'  equations  for  a  thin  shell,  STAGS  uses 
the  Kirchhoff -Love  h)rpothesis  for  strains  away  from  the  midplane. 

With  these  assumptions  in  mind,  an  overview  of  STAGS'  nonlinear 
kinematic  relations  can  now  begin.  The  complete  derivation  of  the 
nonlinear  kinematic  relations  is  given  in  Appendix  A. 

To  start  the  derivation,  the  midsurface  strain,  £°,  is  given  in 
a  form  which  includes  an  in-plane  rotation  term,  <f> ,  [7]: 

-  u,  +  l/2(u,^  +  w,^  +  (2.8) 

X  X  X  X 

In  Sanders'  nonlinear  shell  equations  presented  in  the  last  section 

4>  -  l/2(u,  -  V,  )  (2.9) 

y  X 

which  represent  a  rotation  about  a  normal  to  the  midplane.  In  a 
plate,  the  rotation  of  a  line  segment  within  the  midplane  will  result 
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in  a  vector  normal  to  this  plane,  similar  to  the  rotation  given  in 

Equation  (2.9).  This  normal  rotation  term  for  strain  in  a  flat  plate 

in  the  x-direction  can  be  represented  by  v,  .  This  does  not  imply 

keeping  the  v,^  term  in  Equation  (2.9),  dropping  the  u,  term,  and 

then  substituting  the  results  into  Equation  (2.8).  What  it  says  is 

that  the  normal  rotation  term  is  important  for  use  in  a  plate  in 

order  for  it  to  adequately  model  a  shell.  With  this  in  mind  the 

rotational  term  ^  for  a  flat  plate  is  replaced  directly  by  v,  in 

Equation  (2.8).  Using  this  same  line  of  reasoning,  to  adequately 

track  midsurface  normal  rotations,  the  other  strain  terms,  £°and  -y°  , 

y  xy 

can  be  found.  The  final  expressions  for  the  midplane  strains  are  [7] 


1 

o 

,  +  l/2(u,^  +  V, 

+  w,^) 

X 

X  X 

X  X 

o 

-  V, 

+  l/2(u,^  +  V. 

2  2. 

.  +  W,  ) 

y 

y  y 

y  X 

0 

'  -  U, 

+  V,  +  (u,  u. 

+  V,  V, 

xy  y  X  X  y  X  y 


+  w,  w,  ) 


X  y 


(2.10) 


If  the  Kirchhof f -Love  hypothesis  is  considered  for  out  of  plane 
strain  terms  (see  Appendix  A)  and  combined  with  the  in-plane  strain 
terms  from  Equation  (2.10),  the  resulting  expressions  for  strains  in 
the  plate  are  given  by: 


e 

X 


o 

e 

X 


z  w, 

XX 


6  -  €  -Z  W, 

y  y  yy 


7-7°  -22  w. 


(2.11) 


xy  xy  xy 

Equation  (2.11)  shows  the  nonlinear  kinematic  equations  which 
appear  to  be  used  in  STAGSC-1.  These  kinematic  relations  allow  for 
large  displacements  and  moderate  rotations  (due  to  the  Kirchhof f- Love 
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hypothesis)  and  will  be  used  in  the  following  section  to  derive  a 
general  form  of  the  tangent  stiffness  matrix. 


2.2.2  Derivation  of  the  Tangent  Stiffness  Matrix 

The  tangent  stiffness  matrix  is  a  nonlinear  stiffness  matrix 
used  in  the  Newton-Raphson  method  (see  section  2.2.5)  for  solving  a 
nonlinear  set  of  equations.  There  are  other  techniques  for  solving 
these  equations,  but  STAGSC-1  uses  the  Newton-Raphson  (or  modified 
Newton-Raphson),  therefore,  the  derivation  of  the  tangent  stiffness 
matrix  is  of  interest  in  this  thesis.  The  tangent  stiffness  matrix 
is  derived  for  a  flat  plate  element,  a  STAGS  type  element,  in  a 
general  way  without  reference  to  any  specific  shape  functions.  For 
the  derivation  of  the  specific  shape  functions  used  in  the  STAGS 
element  see  reference  [13].  It  must  also  be  pointed  out  that  STAGS 
uses  an  isoparametric  formulation  whereas  this  formulation  is  given 
in  general  coordinates  to  show  the  reader  the  steps  involved  in 
formulating  a  nonlinear  stiffness  matrix. 

To  start  the  derivation,  a  form  of  the  tangent  stiffness  matrix 
is  found  from  the  energy  expression.  Consider  {'i'),  the  sum  of 
external  and  internal  generalized  forces  (Appendix  B,  Equation 
(B.14)),  which  is  given  as  [8] 

{'I')  -  I  [B]^{c7)  dVol  -  {f}  (2.12) 

•'v 

where 

Iff)  -  vector  of  stresses 

{f)  -  nodal  forces  of  the  element 
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and  [B]  Ls  defined  as  the  product  of  the  differential  operator 
matrix,  [L]  ,  operating  on  the  element  shape  functions,  [N] .  In  an 
equilibrium  state  {^)  (Equation  (2.12))  will  equal  zero.  The  strain 
displacement  relations  can  be  written  (to  include  the  [B]  matrix)  as 
{£)  -  [L][N]{a}  -  [B]{al  (2.13) 

where  (a)  is  the  nodal  displacement  vector.  In  the  case  of  a 
nonlinear  stiffness  matrix  (i.e.  tangent  stiffness  matrix),  [B]  is 
redefined  as  [8] 

[B]  -  [Bo]  +  [Bl]  (2.14) 

where 

[Bo]  -  constant 

[Bl]  -  function  of  displacements 

In  order  to  use  the  Newton -Raphson  method,  a  relation  between  d{a} 
and  dlit)  (see  section  2.2.5)  must  be  found.  Taking  the  variation  of 
Equation  (2.12)  with  respect  do  d{a)  gives  the  relation  needed  and  is 
given  as  [8] 

d{ir)-  I  d[B]'^(a)  dVol  +  |  [B]^d(a)  dVol  -  [KT]d(a)  (2.15) 

J  V  J  V 

where  [Kt]  is  the  tangent  stiffness  matrix.  In  this  derivation, 
strain  is  assumed  small  therefore  the  equation 

(a)  -  [D]  ({£)-  {€  ))  +  {cr  )  (2.16) 

o  o 

still  applies,  where  [D]  is  the  material  matrix  and  the  subscripts 
indicate  initial  values  of  stresses  and  strains  (i.e.  constants). 

With  Equations  (2.13),  (2.14),  and  (2.16),  the  variational  terms  in 
Equation  (2.15)  can  be  rewritten  as  [8] 

d{a)  -  [D]d(£)  -  [D] [B]d{a)  (2.17) 
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and 


d[B]^-  d[BL]^  (2.18) 

With  Equations  (2.17)  and  (2.18),  Equation  (2.15)  can  be  rewritten  as 

dl*)  -  [  d[BL]^cr)  dVol  +|  tB]^[D]  [B]  dVol  d{a)  (2.19) 

jy  Jy 


The  first  term  in  Equation  (2.19)  contains  the  initial  stress 
matrix,  [Ko] ,  while  the  last  term  turns  out  to  be  the  linear  and  the 
nonlinear  stiffness  matrices,  [Ko]  and  [Kl]  respectively,  after 
substituting  Equation  (2.14)  in  for  [B]  [8].  The  three  stiffness 


The  full  expression  for  the  tangent  stiffness  matrix  can  now  be 
written  as 


[Kt]  -  [Ko]  +  [Kl]  +  [Ko]  (2.21) 

With  an  expression  for  the  elements  contained  in  the  tangent 
stiffness  matrix  (Equation  (2.21)),  the  derivation  can  now  proceed 
toward  finding  the  terms  contained  in  the  matrices  that  make  up  the 
tangent  stiffness  matrix. 

The  next  steps  in  the  derivation  include  formulating  the 
kinematic  relations,  introducing  the  material  matrix,  and  giving  the 
displacement  relations.  Recalling  the  strain  expression  from  the 
last  section  (Equation  (2.11))  written  in  vector  form  and  rearranged 
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to  reflect  in-plane  and  bending  (out  of  plane)  strains  [8] 


{€) 


z 

o 

6 

y 

o 

1 

xy 

-w, 

] 

-w, 

-2w, 


yy 


(2.22) 


zy  ' 

where  the  distance  from  the  raidplane,  z,  has  been  incorporated  in  the 
material  matrix  (see  Appendix  C) .  If  Equation  (2.15)  is  rewritten, 
separating  linear  and  nonlinear  terms  due  to  in-plane  and  bending 
strains,  the  results  are  [8] 


{«)  - 


f  •• 

0 

L 

b 

€ 

o 

N.  J 

►  +  • 

0 

>  +  > 

(2.23) 


where 

{e**)  -  linear  in-plane  strains 

O 

{€^)  -  linear  bending  strains 

O 

(e^)  -  nonlinear  in-plane  strains  in  u  and  v 
(6^)  -  nonlinear  in-plane  strains  in  w 

Incorporating  the  linear  constitutive  relationship,  the  material 
matrix  [D] ,  for  a  composite  material  (derived  in  Appendix  C) ,  is 
given  as 


[D] 


[A]  [B]‘ 

[B]  [D] 


(2.24) 


Since  the  distance  from  the  midplane,  z,  has  been  incorporated  in 
this  matrix,  the  integration  indicated  in  Equation  (2.20)  for  [Ko] , 
[Kl] ,  and  [Ka]  reduces  to  an  area  integral. 

Finally,  the  displacements  are  defined  in  terras  of  nodal 
displacements  using  the  shape  functions  for  the  plate  element.  If 
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for  example,  an  element  similar  to  the  QUAF  410  is  considered  (see 
section  2.2.4),  which  has  four  comer  nodes  and  six  degrees  of 
freedom  per  node ,  the  displacements  can  be  given  as 


{ 


[N]{a  1 


(2.25) 


The  vector  of  element  nodal  displacements  (parameters),  (a^),  can  be 
subdivided  into  displacements  that  influence  in-plane  (superscript  p) 
and  bending  (superscript  b)  as  [8] 


(a^)  - 


la^)  - 


u 


W 

i 

w 

L  w,*  j 

yi 


(2.26) 


where  /3  represents  an  average  rotation  about  the  normal  (similar  to  4> 
in  Equation  (2.9)).  The  shape  functions  can  also  be  divided  in  a 
similar  manner  as  [8] 


[N  J 


[N^]  0 

0  [N^]J 


(2.27) 


In  order  to  proceed  further  it  is  necessary  to  expand  the 
expression  for  [B] .  From  Equation  (2.14)  it  was  shown  that  [8] 

[B]  -  [Bo]  +  [Bl]  (2.28) 

where 


[Bo]  . 

(6x24 ) 


and 


[Bl]  - 

(6x24 ) 


[B^]  0 

(3x12) 

0  [b‘iJ 

(3x12) 


[Bj]  [B^] 

(  3x12) ( 3x12) 


(2.29) 


0 


0 


(2.30) 
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Note;  The  numbers  in  the  parenthesis  incate  the  matrix  order. 


The  matrices  [B**]  (linear  in-plane)  and  [B**]  (linear  bending)  are 

O  O 

standard  matrices  and  are  derived  in  reference  [8]  with  the  exception 
that  in  the  derivation  of  [B*"]  the  normal  rotation  term,  ^  (Equation 

O 

(2.26)),  must  be  included  as  a  degree  of  freedom  at  each  node.  The 
matrix  [bI*]  and  [bN  are  found  by  taking  a  variation  of  (c’^)  and  (e*’) 
(Equation  (2.23))  respectively,  with  respect  to  the  nodal  degrees  of 
freedom  (a) . 

In  matrix  form  {€^)  can  be  written  as 


u. 

0  V,  O' 

r 

u 

•v 

_  1 

X 

X 

0 

c 

o 

< 

U 

y 

2 

y  y 

V 

I-, 

U,  V,  V, 

V 

X 

y 

X  y  X  '' 

y  . 

(2.31) 


or  (8] 

{€[)  -  1/2  (R‘’]{S'’)  (2.32) 

The  vector  can  be  related  to  the  in-plane  (u  and  v)  nodal 

degrees  of  freedom  as 


[O' 


u, 

u,; 

V,' 


-  [G^’l  {a®’) 

(  <>xl2  )  (  12x1) 


(2.33) 


where  [G**]  is  a  matrix  of  coordinates  (i.e.  the  derivatives  of  the 

in-plane  shape  functions)  and  (a**)  is  the  vector  of  nodal 

displacements  (Equation  (2.26)).  The  next  step  involves  taking  the 
P 

variation  of  To  do  this  Equation  (2.31)  is  rewritten  as 
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(2.34) 


I 

I 

I 


'•J' 


l/2u,^  +  1/2  v.^ 

l/2u.^  +  1/2  v/ 
y  y 

U,  U,  +  V,  V, 

X  y  X  y 


Taking  the  variation  of  Equation  (2.34)  gives 


d{ef) 


u,  du,  +  V,  dv, 

XX  XX 

u,  du,  +  V,  dv, 
y  y  y  y 

u,  du,  +  u,  du,  +  V,  dv,  +  V,  dv, 
xy  yx  xy  yx 


Rewriting  Equation  (2.35)  results  in 


d{ep  - 


u. 

0 

V, 

0 

u 

X 

X 

X 

U 

0 

u, 

0 

V, 

d  < 

V 

y 

y 

y 

V 

X 

u, 

u, 

V, 

V, 

y  ^ 

y 

X 

y 

X 

(2.35) 


(2.36) 


where  the  first  matrix  on  the  right  hand  side  equals  [R*^]  (Equation 
2.32)  and  the  last  matrix  can  be  written  from  Equation  (2.33)  as 
u, 

:  I  _  _ 

(2.37) 


u, 


d-{  V, 
V, 


[G‘’]d(a‘’) 


( 4x  12  )  (  12x  1) 


since  [G**]  is  a  function  of  coordinates  only. 

Using  the  definition  of  [R**]  in  Equation  (2.32)  along  with 
Equation  (2.37),  Equation  (2.36)  can  be  written  as 


d{£^)  -  [R*’]  [G'’]d{a‘’) 

From  Equation  (2.38)  the  resulting  expression  for  [B^]  is,  by 
definition, 


(2.38) 


[B^]  -  (R‘']  [G**] 

(3x12)  ( 3x4 ) (4x12) 


(2.39) 
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Following  the  same  line  of  reasoning  used  to  find  [B^]  from 
one  can  find  from 


(2.40)  ! 

I 


or  from  [8] 

-  1/2  (2.41) 

that 

[b"]  -  [r"]  [g"]  (2.42) 

(3x12)  (3x2)  (2x12) 

where  [G^]  is  composed  of  derivatives  of  the  shape  functions  that  are 
contained  in  the  expression  for  w. 

With  the  expressions  for  [Bo]  and  [Bl]  determined  (Equations 
(2.29)  and  (2.30)  respectively),  the  material  matrix  defined 
(Equation  (2.24)),  and  recalling  that  the  volume  integral  has  been 
reduced  to  an  area  integral  (Equation  (2.24)),  the  expressions  for 
the  linear  and  nonlinear  stiffness  matrices  ([Ko]  and  [Kl])  can  be 
determined.  If  one  uses  Equation  (2.20)  and,  after  some 
manipulation. 


[Ko] 

(2<.  xZit  ) 


■I 


[b‘’]^[A][B‘’]  [b‘’]^[B](b‘’] 

o  o  o  o 

[b'’]^[B][b‘’]  [b‘’]^[D][b’’] 

o  o  o  o 


dArea 


and 


[Kl] 

(24x24 ) 


■I 


[1] 

[2] 

[2] 

[3] 

dArea 


(2.43) 


(2.44) 


where 
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(2.45) 


[1]  -  [B‘’]^[A][Bf]  +  [Bf]^[Al[Bf]  +  [Bf]^[A][B'’] 

O  L  L  L  L  O 

(12x12) 


[2]  -  [bP]"[A][B;’]  +  [Bf]"[A]{B^  +  {bJ]"[B][B^ 

O  L  L  L  L  O 

(12x12) 


(2.46) 


[3]  -  [b‘’]"[B][B;’]  +  [bJ’i'IAIEB^]  +  [Bh"[B][B^ 

O  L  L  L  L  0 

(12x12) 


(2.47) 


The  final  expression  necessary  for  the  tangent  stiffness  matrix 
is  the  initial  stress  matrix  [Ka],  Recalling  Equation  (2.20), 


rewritten  for  convenience  as 


[Ko]d(a) 


I  <i[BL]^^  ' 


(2.48) 


The  stresses,  [a  '],  are  defined  in  terms  of  the  in-plane 
(superscript  p)  and  bending  (superscript  b)  stresses  as  [8] 


[o')  - In  n  nm  m  m1  • 

i  X  y  xy  X  y  xyj 


(2.49) 


where  the  prime  on  the  stresses  indicates  stress  resultants,  not  true 
stresses,  and  the  true  stress  resultants  (N  ,  M  ,  etc.)  are  given  in 

XX 

Appendix  C  (Equation  (C.ll)).  The  stress  resultants  are  average 
values  over  the  element.  This  also  implies  that  integration  through 
the  thickness  has  been  completed  thereby  reducing  the  integration  of 
Equation  (2.48)  to  that  of  an  area.  Now,  define  the  variation  of 
[Bl]^  in  Equation  (2.48)  from  Equation  (2.30)  as 


dlB^]"  0 
d[B^]"  0 


(2.50) 


Substituting  Equations  (2.39)  and  (2.42)  into  Equation  (2.50)  and 
then  substituting  that  result,  as  well  as  Equation  (2.49),  into 
Equation  (2.48)  gives  [8] 
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[Ka]d{a}  -J 

A 


■[G*’]^d[R*’]*  O' 
_[G‘’]*d[R‘’]^  0_ 


(2.51) 


Expanding  Equation  (2.51)  and  rewriting  it  in  terms  of  in-plane 
(superscript  p)  and  bending  (superscript  b)  expressions  results  in 


[Ka*’]d{a‘’} 


.[  [G^]^  d[R'’]^ 


(12x 12 )  ( 12x  1  )  "(12x4)(4x3) 


dArea 


(2.52) 


[Ka‘’]d{a‘’)  -J  [G‘’]^d[R'']^-  N  •  dArea 

(12xl2)(12xl)  •*(12x2X2x3) 


(2.53) 


The  following  steps  involve  finding  expressions  for[Ka*’]  and  [Ka^]  . 

Starting  with  Equation  (2.52)  and  rewriting  the  last  two 
matrices  on  the  right  hand  side  using  Equation  (2.31)  gives 

u,  0  u, 

X  y 

0  u,  u, 

y  X 

V,  0  V. 

X  y 

0  V,  V, 

y  X 

Expanding  the  right  hand  side  of  Equation  (2.54)  and  also  taking  the 
variation,  recalling  that  N  ,  N  ,  and  N  are  constants  (Equation 

X  y  xy 

(2.49)),  results  in 
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N 

N  Y  - 
y 

xy^ 


f  N  du,  +  N  du,  1 

XX  3cy  y 

N  du,  +  N  du, 
y  y  xy  X 

N  dv,  +  N  dv, 

XX  Tf  y 

N  dv,  +  N  dv, 
y  y  xy  X 


(2.55) 


Rewriting  Equation  (2.55)  in  matrix  form  gives 


d[R 


PiT 


■  N 

N 

0 

0 

r  N  1 

X 

X 

N 

xy 

N 

0 

0 

u, 

u. 

•V 

'  X 

>  N 

xy 

y 

d< 

V, 

y 

y 

0 

0 

N 

N 

V, 

X 

N 

xy^ 

X 

xy 

1. 

y  - 

0 

0 

N 

N 

. 

xy 

y 

(2.56) 


Recalling  Equation  (2.37)  and  substituting  it  into  Equation  (2.56) 
gives 


d[R' 


PiT 


N 

X 

N,h- 

”xy^ 


N  N  0  0 

X  xy 

N  N  0  0 

xy  y 

0  0  N  N 

X  xy 

0  0  N  N 

xy  y 


[G^]d{a^) 


(2.57) 


Substituting  Equation  (2.57)  into  Equation  (2.52)  results  in  an 
expression  for  [Ka’’]: 


[Ka*’]  - 

(12x12) 


[G'^r 

(12x4) 


N  N  0  0 

X  xy 

N  N  0  0 

xy  y 

0  0  N  N 

X  xy 

0  0  N  N 

xy  y 


[G  ]  dArea 
(4x12  ) 


(2.58) 


With  the  same  type  of  derivation  used  to  find  [K<7^]  ,  [Ka^]  can 


be  found.  The  resulcing  expression  for  [Kcr  ]  is: 
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■  N 

N  1 

[Ka^]  - 

[G*”]’ 

X 

xy 

[G^]  dArea 

(2.59) 

(12x12)  , 

(12x2) 

N 

N 

(2x12 ) 

i 

i 

^  xy 

y 

With  Equations  (2.58)  and  (2.59),  the  full  expression  for  the 
initial  stress  matrix  can  be  written  as 


[Ka]  - 

(2Ax24  ) 


[KO**]  0 

(12x12) 

0  [Ko*"] 

(12x12)-* 


(2.60) 


2.2.3  Pressure  Loading 

Since  STAGSC-1  is  an  energy  based  finite  element  program,  the 
loading  is  restricted  to  conservative  systems.  In  a  conservative 
force  system  the  work  done  by  these  forces  is  path  independent  [13] . 
In  this  thesis  the  loading  is  an  internal  pressure  load, 
asymmetrically  applied.  A  pressure  load  like  this  that  remains 
normal  to  the  surface  during  deformation  (a  live  load)  is  not 
necessarily  conservative  [13].  An  example  of  a  conservative  pressure 
system  is  a  pressure  acting  on  a  closed  body;  the  work  equals  the 
product  of  the  pressure  and  the  change  in  volume.  For  plates  and 
shells,  the  pressure  loading  is  conservative  if,  at  the  edge  of  the 
pressure  field,  the  normal  displacement  or  the  product  of  the 
displacement  vector  (u,v)  and  the  normal  to  the  boundary  in  the 
tangent  plane  is  zero  [7,15]. 

From  extensive  research  on  STAGSC-1  theory  and  other  related 
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articles  on  pressure  loading,  it  would  appear  that  the  live  pressure 
loading  in  STAGS  is  implemented  as  follows.  The  local  reference 
frame  of  the  element  is  established  through  the  use  of  the  updated 
Lagrangian  formulation  (see  section  2.2.5).  A  normal  is  calculated 
to  this  local  element  reference  frame  and  is  established  as  the 
direction  by  which  the  pressure  load  is  applied.  Equivalent  nodal 
loads  are  calculated  for  the  element  for  the  current  pressure  load 
(loads  are  incrementally  applied).  These  loads  are  transformed  into 
the  global  system  and  assembled  with  the  other  elements'  equivalent 
nodal  loads.  A  solution  is  then  found  for  the  displacements  at  this 
given  load  increment  (the  global  force  vector  doesn't  change  within  a 
load  increment) .  Once  these  displacements  are  found  ,  the  element 
reference  system  as  well  as  the  element  normal  can  be  updated,  an 
increment  of  pressure  load  added,  and  the  process  repeated. 

2.2.4  Plate  Elements  Representing  a  Shell 

As  stated  earlier,  STAGS  uses  plate  elements  to  represent  the 
surface  of  a  shell.  The  quadrilateral  plate  elements  primarily  used 
in  this  thesis  are  the  STAGS  QUAF  410  and  411  elements.  The  QUAF  410 
element  has  24  degrees  of  freedom;  three  rotations  and  three 
displacements  at  each  corner  node.  The  QUAF  411  element  has  32 
degrees  of  freedom;  four  rotations  and  three  displacements  at  each 
corner  node  and  an  in-plane  tangent  displacement  at  each  of  four  side 
nodes .  These  Elements  are  shown  in  Figure  2.2. 

From  examining  Figure  2.2  it  is  noted  that  a  degree  of  freedom 
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not  usually  used  in  plate  elements  is  included;  the  normal  rotation 


d  .  This  normal  rotation,  for  the  QUAF  410  and  411  elements,  is  the 

z 

average  in-plane  rotation  of  the  two  adjacent  edges  of  the  plate 
element.  This  degree  of  freedom  is  necessary  when  two  flat  plate 


/e,.  /e^  /e„  /e,. 

Figure  2.2  STAGS  QUAF  410  and  411  Quadrilateral  Element  [16] 

elements  meet  at  an  angle  to  represent  a  curved  shell  [13],  A 
problem  with  this  rotation  is  that  it  does  not  appear  in  the  strain 
energy  expression  for  the  element. 

To  illustrate,  consider  two  flat  plate  elements  joined  together 
at  an  angle  a  as  in  Figure  2.3.  These  plates  represent  the  curved 
surface  of  a  shell.  Once  the  rotations  of  each  element  are 
transformed  into  the  same  reference  system,  compatibility  is  given  as 
[13] 

jcos(Q/2)  -  [  jsin(Q/2)  -  0 


where  the  superscripts  1  and  2  are  the  associated  element  numbers. 

As  the  angle  between  the  elements,  a,  becomes  smaller  the  system  of 
equations  in  Equation  (2.61)  becomes  increasingly  ill-conditioned  and 


Figure  2.3  Flat  Plate  Element  Representing  a  Curved  Shell  [13] 

is  singular  at  a-0  [13],  In  order  to  overcome  this  difficulty  a 

lower  limit  is  set  for  the  angle  a.  Once  this  limit  is  reached  9  is 

I  z 

omitted  and  Equation  (2.61)  is  assumed  as  ^  . 

y  y 

Another  problem  associated  with  the  use  of  flat  elements 
representing  a  curved  surface  has  to  do  with  displacement  conform!  ty 
at  an  interface  of  two  flat  elements.  For  a  flat  element  the  lateral 
deflection,  w,  is  usually  represented  by  at  least  a  cubic  polynomial 
in  order  to  handle  the  second  derivatives  associated  with  bending 
strain.  The  in-plane  displacements,  u  and  v,  are  usually  represented 
with  either  a  linear  or  quadratic  polynomial.  If  the  compatibility 
relation  for  two  adjacent  flat  elements  is  derived,  the  resulting 
expression  is  given  as  [13] 

f  v^  -  v^  ]cos(a/2)  -  f  w^  +  w*  ]sin(a/2)  -  0  (2.62) 
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jcos(a/2)  -  j  jsin(a/2)  -  0  (2.62) 

^  ^  ^  (cont) 

If  one  examines  Equation  (2.62),  it  is  evident  that  the  displacement 

compatibility  along  the  interface  of  the  two  elements  cannot  be 

satisfied  (if  v  is  a  quadratic,  it  cannot  fit  the  same  curve  as  the 

cubic,  w)  unless  v  and  w  are  represented  by  polynomials  of  the  same 

order.  To  overcome  this  displacement  nonconformity  the  QUAF  410  uses 

a  third  order  polynomial  to  represent  u  in  the  y-direction  (linear  in 

x)  and  V  in  x-direction  (linear  in  y) .  The  lateral  deflection  w  is 

represented  by  a  cubic  polynomial  in  x  and  y.  There  are  a  total  of 

twelve  terms  in  polynomials  representing  u  and  v  for  the  QUAF  410 

element.  This  can  be  done  since  there  are  two  displacements,  u  and 

v,  and  one  rotation,  9  ,  at  each  node  representing  in-plane  motion. 

2 

The  QUAF  411  element  uses  the  same  cubic  polynomial  to  represent  w 
but  adds  a  shear  term,  7,  at  each  corner  node  and  an  in-plane  tangent 
displacement,  t,  at  midside  nodes.  This  allows  u  to  be  cubic  in  y 
(quadratic  in  x)  and  v  to  be  cubic  in  x  (quadratic  in  y)  [13]. 
Rotational  compatibility  is  enforced  only  at  the  nodes,  thus  the  QUAF 
410  and  411  elements  are  nonconforming  bending  elements.  The 
derivation  of  the  shape  functions  associated  with  these  assumed 
displacement  fields  is  given  in  [13]. 


2.2.5  Solution  Techniques 

The  STAGSC-1  computer  program  uses  a  Newton-Raphson  or  modified 
Newton- Raphson  (user  defined)  solution  technique  to  solve  the 


26 


nonlinear  equilibrium  equations.  For  a  linear  problem  in  STAGSC-1, 
the  solution  technique  is  much  easier  since  the  stiffness  matrix  is 
not  a  function  of  displacements.  The  linear  set  of  equations  is 
given  as 

[K]{a)-{R}  (2.63) 

where  [K]  is  the  stiffness  matrix  (constant),  (a)  is  the  set  of  nodal 
degrees  of  freedom  (nodal  parameters),  and  (R)  is  the  applied  loads. 
STAGSC-1  uses  a  Cholesky  triangular  decomposition  with  forward  and 
backward  substitution  to  solve  Equation  (2.63)  for  (a)  [12].  The 
solution  for  the  nonlinear  equilibrium  equations  is  much  more 
involved  and  will  be  looked  at  next. 

Before  starting  the  explanation  of  the  Newton-Raphson  solution 
technique,  an  explanation  of  the  reference  systems  in  STAGSC-1  is  in 
order.  There  are  basically  two  reference  systems  defined;  a  global 
reference  system  (subscript  g)  and  a  local  element  reference  system 
(subscript  e) .  The  global  reference  system  is  fixed  in  space  and 
does  not  move.  The  local  reference  system  uses  an  updated  Lagrangian 
approach.  The  local  reference  system,  called  a  corotational  system, 
is  fixed  to  the  element  and  moves  with  the  element  during  rigid  body 
motion  [6].  This  allows  for  the  removal  of  rigid  body  motion  in  the 
element  before  calculating  strain.  The  present  version  of  STAGSC-1 
(1986)  also  redefines  the  standard  way  of  representing  a  rotation  as 
a  vector  quantity.  Rotations  are  defined  by  a  triad  (three  mutually 
perpendicular  unit  vectors)  to  accurately  map  local  rotations  since 
they  are  dependent  on  order  of  rotation.  The  previous  versions  of 
STAGS  used  vectors  to  describe  rotations  which  with  large  rotations 
gave  spurious  results  (vectors  cannot  track  order  of  rotation  and  are 
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only  good  for  small  rotations) .  This  is  the  reason  the  QUAF  411 
element  was  developed;  to  handle  larger  rotations,  but  with  increased 
cost  [18].  The  reader  is  referred  to  [17]  for  a  complete  discussion 
on  this  rotational  formulation.  Finally,  differentiation  and 
integration  within  the  element  are  done  with  respect  to  this 
corotational  system  [6] . 

For  the  STAGSC-1  quadrilateral  element,  the  local  reference 
system  is  defined  as  follows:  An  approximation  of  the  warped  element 
surface  is  made  by  crossing  the  principal  diagonals  of  the  element  to 
form  a  vector  and  then  establishing  a  plane  normal  to  this  vector 
such  that  one  of  the  nodes  lies  on  this  plane.  This  node  is  where 
the  local  reference  system  is  established  by  taking  the  z  axis  normal 
to  the  plane,  the  x  or  y  axis  is  located  along  one  of  the  element 
edges,  and  the  remaining  axis  completes  a  cartesian  right-handed 
system  (see  Figure  2.4)  [17].  All  local  deformations  of  the  element 
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1 


Figure  2.4  Quadrilateral  Element  Local  Reference  System  [17] 

are  referenced  to  this  local  system  during  an  iteration  for  a 
solution.  Once  a  solution  has  been  found  in  the  global  system  the 
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local  coordinates  are  updated  and  a  new  local  reference  system  is 
established.  This  movement  of  the  local  reference  system  appears  to 
be  Eulerian  except  that  the  local  coordinates  of  a  point  change  [6] . 
Strains  and  rotations  in  the  local  reference  system  are  usually  small 
enough  so  as  to  ignore  the  nonlinear  stiffness  matrix  [Kl]  and 
sometimes  even  the  initial  stress  matrix  [Ko]  that  are  part  of  the 
tangent  stiffness  matrix  [Kt]  [6] .  This  is  only  true  for  the  local 
system,  the  global  tangent  stiffness  matrix  must,  in  general,  contain 
all  terms  (see  Equation  (2.21)).  From  thorough  research  of  STAGS 
documentation,  it  appears  that  STAGS  uses  the  full  tangent  stiffness 
matrix  in  both  the  global  and  local  systems  to  overcome  difficulties 
in  solutions  due  to  highly  nonlinear  problems. 

With  the  reference  systems  defined,  the  Newton- Raphson  solution 
techniques  can  be  explained.  To  start,  expand  the  expression  for  the 
sum  of  external  and  internal  forces,  {'P)  (see  Appendix  B) ,  in  a 
Taylor  expansion 

|'I-({a)"*^)j  -j^f({a)")J-  +  (Aa)"  +H.O.T.  -  0  (2.64) 

'■  ■'  n 

where  H.O.T.  means  higher  order  terms  and 

(a)"''^  -  (a)"  +  (Aa)"  (2.65) 

From  Equation  (2.15)  [8] 

(2.66) 

where  [Kt]  is  the  tangent  stiffness  matrix;  a  function  of 
displacements.  If  one  uses  the  expression  for  {4'}  in  Equation 
(B.15)  (Appendix  B)  and  writes  it  as  a  function  of  displacements  for 
the  nonlinear  problem  one  obtains  [8] 


{«P(a))  -  [K(a)]{a)  -  (f)  (2.67) 

where 

[K(a)]{a)  -  internal  resisting  forces  of  the  element 
If)  -  externally  applied  loads 

Inserting  Equations  (2.66)  and  (2.67)  into  Equation  (2.64)  and 
eliminating  the  higher  order  terms  results  in 

[KT]{Aa}"-  If)  -  [K(a)"](a)"  (2.68) 

The  problem  now  is  to  find  the  the  displacements  within  the  element 
that  balance  the  externally  applied  loads  and  the  internal  resisting 
forces  (i.e.  the  sum  of  the  right  hand  side  of  Equation  (2.68) 
equaling  zero) .  Since  an  updated  Lagrangian  approach  is  being  used, 
the  internal  deformations  are  relative  to  the  local  reference  system, 
therefore  the  last  expression  in  Equation  (2.68)  needs  to  be  written 
in  terms  of  the  local  displacements.  Rewriting  Equation  (2.68)  in 
terms  of  the  local  (subscript  e)  and  global  (subscript  g) 
displacements  results  in  [6] 

[KTllAa^)"  -  (f)  -  X  [k(a^)"]{aj"  (2.69) 

where 

^  [k(a  )"]{a  )"  -  internal  resisting  forces  transformed 
into  the  global  system  and  summed 

The  Newton-Raphson  solution  technique  is  given  as  follows  [6]: 

1.  Increment  the  applied  load  {f). 

2.  Establish  the  local  coordinates  for  each  element. 

3.  Formulate  the  element  tangent  stiffness  matrices  in  terms  of 
their  local  degrees  of  freedom. 

4.  Transform  the  local  tangent  stiffness  matrices  to  the  global 
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coordinate  system  and  assemble  them  Into  the  global  stiffness 
matrix. 

5.  Compute  the  values  of  the  local  degrees  of  freedom  (a^) 

(zero  for  the  first  Iteration  within  each  load  step)  from  the 
global  degrees  of  freedom 

6.  Calculate  the  element  Internal  forces  [k(a  )]{a  ). 

•  9 

Transform  these  forces  to  the  global  system  and  assemble  them 
with  the  other  element  Internal  forces. 

7.  Solve  Equation  (2.69)  for  {da  ). 

8.  If  the  vector  {da  )  Is  not  small  enough  (l.e.  converges) 
return  to  step  3. 

After  the  solution  converges  (step  8)  return  to  step  1  and  repeat. 
Once  the  solution  has  converged  for  the  final  load  step  ({f)  has  been 
Incremented  to  equal  the  final  load)  the  solution  Is  complete. 

An  alternate  solution  technique  Is  the  modified  Newton-Raphson 
method.  In  this  method  the  tangent  stiffness  matrix  Is  not 
reformulated  at  every  Iteration  within  a  load  step.  The  previously 
assembled  stiffness  matrix  Is  used  for  successive  iterations  and  Is 
only  refactored  when  the  convergence  rate  dictates  [6] . 

STAGS  allows  the  user  to  control  many  of  the  parameters  in  the 
solution  process  to  include:  either  using  the  full  or  modified 
Newton-Raphson  method,  control  of  the  load  step  size,  and  the  number 
of  attempts  at  a  solution  if  convergence  difficulties  arise. 
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3.  Finite  Element  Modeling 

The  analysis  carried  out  in  this  thesis  consists  of  a  linear  and 
nonlinear  static  analysis  using  the  Structural  Analysis  for  General 
Shells,  version  C-1,  (STAGSC-1)  finite  element  computer  program. 

The  structure  modeled  is  a  thin  composite  shell  with  the 
following  characteristics: 

1.  Varying  thickness  over  the  shell. 

2.  Kevlar-49/F-141  Polyester  Resin  (Kevlar/Polyester)  cloth 

composite . 

3.  Dimensions  of  approximately  216  in.  long,  38  in.  wide,  and 

34  in.  high. 

4.  Unusual  boundary  conditions  (not  truly  clamped  or  pinned). 

5.  Internal  asymmetric  pressure  loading  (see  section  2.2.3). 
These  characteristics  will  be  expanded  upon  in  the  following 
sections.  The  actual  shell  is  shown  in  Figure  3.1.  The  primary  area 
of  concern  in  this  analysis  is  the  center  section  of  the  shell  away 
from  the  boundaries;  the  area  where  experimental  results  were 
obtained  (see  section  4). 

To  model  the  structure,  a  finite  element  grid  was  developed  to 
approximate  the  curved  surface  of  the  shell  with  flat  quadrilateral 
and  triangular  plate  elements.  The  basic  model  consists  of  362  nodes 
and  362  elements.  In  areas  where  geometry  dictates,  smaller  elements 
were  used  to  better  approximate  the  cuirvatures.  The  finite  element 
model  of  the  structure  is  shown  in  Figures  3.2,  3.3,  and  3.4. 

In  order  to  model  this  kind  of  structure  using  STAGSC-1, 
individual  node  point  locations  had  to  be  entered  in  the  global 
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Figure  3.1  Kevlar/Polyester  Composite  Shell 


system  (reference  located  at  the  nose  of  the  model).  Also,  element 
connectivity  had  to  be  entered  by  hand.  This  type  of  model  is  called 
an  element  unit  in  STAGS  [16].  If  the  model  was  a  standard  geometric 
shape  (e.g.  cylinder,  ellipsoid,  etc.),  STAGS  would  generate  the 
finite  element  mesh. 


Figur 
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Figure  3.3  Finite  Element  Model  Front  and  Rear  View 


Figure  3.4  Finite  Element  Model  Angle 
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automatically  and  reduce  the  number  of  user  inputs  dramatically. 

This  standard  geometric  shape  in  STAGS  is  called  a  shell  unit  [16], 
The  loading  on  the  structure  is  an  internal  asymmetric  pressure 
load  over  the  surface  of  the  shell.  The  loads  are  due  to  aerodynamic 
pressures  on  the  shell  and  static  equivalents  were  calculated  by  the 
4950  Test  Wing,  Wright-Patterson  AFB,  Ohio.  The  static  loads  were 
divided  into  20  loading  regions  on  the  shell;  symmetric  with  respect 
to  the  x,y  plane  (see  Figures  3.2,  3.3,  and  3.4)  in  terms  of 
location,  but  not  in  terms  of  load  applied.  Three  different  load 
conditions  were  calculated;  the  "worst  case  loading"  is  presented  in 
this  thesis.  The  loading  regions  are  shown  in  Figure  3.5, 
superimposed  on  the  finite  element  grid.  The  loading  associated  with 
these  regions  is  presented  in  Table  3.1.  This  type  of  loading  was 
achieved  in  STAGSC-1  by  the  use  of  a  subroutine  called  UPRESS  that  is 
user  generated.  The  user  generates  the  code  along  the  guidelines  in 
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Table  3 . 1  Static  Equivalent  Loads 


reference  [16],  compiles  it,  and  links  it  with  the  STAGS  program. 
Since  the  loading  in  this  case  is  a  live  (remains  normal  to  the 
surface)  pressure  load,  the  subroutine  allows  for  a  flag  to  the  main 
program  to  calculate  the  loads  as  follower  (live)  loads  (see  section 
2.2.3).  The  subroutine  used  on  the  final  model  is  shown  in 
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Element  Model  and  Loading  Regions 


Appendix  H.  The  loads  on  the  structure  are  applied  in  0.1  increments 
of  the  total  load  for  nonlinear  runs. 

The  modeling  of  this  shell  structure  was  done  in  three  phases . 

In  each  of  these  phases  the  previous  model  was  modified  to  more 
closely  approximate  the  actual  shell.  For  example,  changes  in  the 
first  phase  were  kept  for  Phase  II  and  improved  upon.  Some  of  the 
modeling  parameters  that  were  changed  to  more  accurately  model  the 
shell  and  the  phase  in  which  they  were  changed  include: 

1.  Thickness  variation  (Phase  I). 

2.  Material  modeling  (Phase  II  and  III). 

3.  Boundary  conditions  (Phase  II). 

4.  Element  type  (Phase  III). 

The  baseline  model  (starting  model)  consists  of  a  constant 
thickness,  isotropic  shell  with  clamped  boundary  conditions  and 
STAGS'  QUAF  410  quadrilateral  and  TRINC  320  triangular  elements  (see 
section  2.2.4).  The  primary  type  of  analysis  during  these  phases  is 
a  static  nonlinear  analysis  with  linear  static  runs  completed  for 
comparison  purposes.  For  the  nonlinear  runs  a  full  Newton-Raphson 
solution  technique  (see  section  2.2.5)  was  used  to  obtain  the 
displacement  solutions. 

Load  versus  normal  (to  the  shell  surface)  displacement  plots  are 
shown  for  Phase  I  through  III  models  at  nodes  163  and  176  (see  Figure 
3.5).  In  the  element  unit  model  the  degrees  of  freedom  coincide  with 
global  axis  directions.  In  order  to  later  compare  the  experimental 
displacements  (normals)  to  the  finite  element  displacements  the 
finite  element  model  displacements  needed  to  be  transformed  into 
surface  normals.  The  method  used  to  transform  the  finite  element 
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global  displacements  to  surface  normal  displacements  is  shown  in 
Appendix  D. 

To  insure  that  the  finite  element  model  accurately  models  the 
shell,  a  convergence  test  was  conducted  to  insure  that  the  aspect 
ratios  used  in  the  model  are  acceptable.  The  aspect  ratio  is  the 
length  of  the  longest  side  of  the  element  divided  by  its  short  side. 
In  general,  elongated  (shaped  like  a  rectangle)  linear  elements 
behave  poorly,  but  quadratic  elements  are  well  behaved  when  elongated 
[6],  Since  the  STAGS  QUAF  410  element  is  cubic  in  one  direction  and 
linear  in  the  other  (see  section  2.2.4),  it  is  not  apparent  whether 
the  QUAF  410  element  will  behave  correctly  in  this  model. 

In  the  center  section  of  the  model,  away  from  the  boundaries, 
the  aspect  ratios  are  approximately  1.5  to  1.  A  convergence  test  was 
done  on  a  similar  structure  with  aspect  ratios  of  2  to  1  and  1  to  1 
to  determine  if  the  model  is  acceptable.  The  following  section 
contains  the  results  of  this  test. 

The  following  sections  detail  the  convergence  test  and  Phases 
I-III  modeling.  With  the  basics  of  the  model  discussed,  the 
explanation  of  the  finite  element  modeling  can  proceed. 


3.1 


Test 


The  convergence  test  for  the  finite  element  model  is  done  on  a 
structure  similar  to  the  actual  model.  The  reason  the  analysis  was 
done  on  a  similar  model  rather  than  the  actual  model  is  that  in  order 
to  reduce  the  grid  size  for  the  actual  model,  new  nodes,  and  element 
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connectivity  would  have  to  be  defined  and  entered  by  hand  (the  basic 
model  consists  of  approximately  750  lines  of  input  while  the  final 
model  has  over  2000  lines).  Using  a  STAGSC-1  ellipsoid  shell  unit 
(see  reference  [16])  on  the  other  hand  reduces  the  input  dramatically 
(approximately  70  lines)  while  still  being  able  to  show  the  validity 
of  the  aspect  ratio  used  in  the  real  model.  Convergence  tests  using 
2  to  1  and  1  to  1  aspect  ratios  are  used  to  validate  the  1.5  to  1 
aspect  ratio  used  in  the  actual  model. 

Before  showing  the  convergence  test,  an  explanation  of  the  STAGS 
shell  unit  reference  system  needs  to  be  addressed.  The  reference 
system  is  based  on  a  surface  system  where  the  shell  is  defined  by 
surface  coordinates  X  and  Y.  The  global  system  coordinates  are  x,  y 
and  2.  Figure  3.6  shows  the  relation  of  the  surface  and  global 
reference  systems  for  the  STAGSC-1  ellipsoid  shell. 


X,Y  measured  in  deg. 


Figure  3.6  Ellipsoid  Surface  (X,Y)  and  Global  (x,y,z) 
Reference  Systems 
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For  the  convergence  test  Rx  is  108.0  in.  and  Ryz  is  26.0  in. 

(see  Figure  3.6)  which  will  approximate  the  length  of  the  shell  but 
will  make  the  convergence  model  wider  than  the  actual  model. 

The  thickness  of  the  convergence  model  is  taken  to  be  constant 
over  the  entire  shell  and  is  0.1455  in.  This  is  the  thickness  used  in 
the  Phase  I  constant  thickness  model  (see  section  3.2). 

Although  the  actual  model  is  made  of  a  Kevlar/Polyester 
composite  fabric,  the  convergence  model  will  treat  the  fabric  as  an 
orthotropic  laminate  with  ply  orientations  of  0  and  90  degrees  . 

There  are  a  total  of  40  plies  used  in  this  model  that  make  up  the 
thickness  of  the  shell.  The  material  properties  for  this  material 
are ; 

-  4.89x10^  psi 

-  4. 23x10^  psi 

i/  -  0.44 
12 

1/  -  0.038 
21 

G  -  1.9x10^  psi 

This  type  of  modeling  treats  the  material  orthotropically  as  in  Phase 
II  Modeling  (section  3.3),  and  is  explained  in  that  section. 

The  finite  element  models  for  the  convergence  test  are  shown  in 
Figure  3.7  with  the  loading  regions  superimposed  on  the  grids.  The 
loads  used  are  the  same  as  those  shown  in  Table  3.1.  The  application 
areas  of  the  loads  varies  from  the  actual  model  but  is  sufficiently 
close.  Once  again,  the  thrust  of  this  convergence  test  is  to 
validate  the  aspect  ratio  used  in  the  actual  model.  The  boundary 
conditions  for  the  model  are  assumed  clamped  at  the  bottom  edge. 

The  convergence  models  were  ran  using  the  nonlinear  solution 
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Figure  3.7  Convergence  Test  Finite  Element  Models  with 
Loading  Regions 


option  in  STAGSC-1.  The  element  used  is  a  QUAF  410  (see  section 
2.2.4)  STAGS  element.  The  data  deck  used  on  the  1  to  1  convergence 
model  is  shown  in  Appendix  E.  The  2  to  1  model  input  deck  is 
identical  to  the  1  to  1  model  deck  except  for  the  N4  card  where  every 
other  X  valve  is  omitted  (see  Appendix  E) .  Loading  on  the  structure 
is  achieved  by  using  the  UPRESS  subroutine  option  in  STAGSC-1  [16]. 
The  user  creates  the  subroutine,  compiles  it,  and  then  links  it  with 
the  main  program.  The  subroutine  is  shown  in  Appendix  F. 

The  convergence  test  uses  the  stresses  in  the  elements  for  model 
comparisons.  Stresses  were  chosen  since  the  strains  are  calculated 
from  displacements,  and  stresses  from  the  strains.  If  there  is  any 
discrepancies  between  models,  they  will  be  magnified  in  the  stress 
calculations.  Stresses  in  the  X  and  Y  directions  (surface 
coordinates)  were  calculated  at  the  element  centroids  for  the 
innermost  composite  ply  using  the  STAGSC-1  post  processor  (POSTP) 
[16]. 

To  compare  the  1  to  1  and  2  to  1  aspect  ratio  models,  stresses 
in  the  X  and  Y  directions  (surface  reference)  were  used  at  X  equals 
75  and  105  degrees.  For  the  2  to  1  aspect  ratio  model  this  X 
coordinate  coincides  with  the  centroid  of  the  elements  shown  in 
Figure  3.7;  the  location  where  the  stresses  were  calculated.  For  the 
1  to  1  aspect  ratio  model  the  stresses  at  the  centroids  of  the 
elements  adjacent  to  X  equals  75  and  105  degrees  were  averaged  and 
used  for  comparison  to  the  2  to  1  aspect  ratio  model. 

The  results  of  this  comparison  for  the  stress  in  the  X  direction 
are  shown  in  Figure  3.8  for  X  equals  75  degrees  and  in  Figure  3.9  for 
X  equals  105  degrees.  The  stresses  in  the  Y  direction  are  shown  in 
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Figure  3.10  for  X  equals  75  degrees  and  in  Figure  3.11  for  X  equals 
105  degrees.  The  results  in  Figures  3.9  to  3.11  are  shown  for  all 
elements  in  the  Y  direction  (measured  in  degrees  circumfrentially) 
for  the  respective  X  value. 

From  examining  Figures  3.8  to  3.11  it  is  evident  that  the  aspect 
ratio  of  1.5  to  1  used  in  the  actual  model  should  approximate  the 
shell  surface  correctly.  This  conclusion  is  based  on  the  fact  that 
the  stress  results  for  the  1  to  1  and  2  to  1  models  are  extremely 
close  in  the  both  the  X  and  Y  directions  over  the  entire  area 
examined. 
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Figure  3.8  Convergence  Test  Results  for  o  at  X-IS  degrees 
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Figure  3.9  Convergence  Test  Results  for  a  at  X-105  degrees 
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Figure  3.10  Convergence  Test  Results  for  a  at  X-75  degrees 
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Figure  3.11  Convergence  Test  Results  for  o^at  X-105  degrees 


3 . 2  Phase  1  Modeling 


Phase  I  consists  of  three  models,  developed  sequentially,  to 
better  understand  the  STAGSC-1  code  and  analysis  involved.  Table  3.2 
presents  the  models  developed  in  Phase  I. 


Mod 

Solution 

Thickness 

Material 

B.C.s 

DOF(ADOF)* 

Elem 

1 

Linear 

Constant 

Isotrop . 

Clamp 

410 

2888(1908) 

362 

2 

Nonlin. 

Constant 

Isotrop . 

Clamp 

410 

2888(1908) 

362 

3 

Nonlin. 

Variable 

Isotrop. 

Clamp 

410 

2888(1908) 

362 

*  DOF”  Degrees  of  Freedom  Elem”  Total  #  of  Elements 


ADOF”  Active  Degrees  of  Freedom  Type”  Element  Type 

Table  3.2  Phase  I  Models 


The  starting  model  in  the  analysis  is  a  constant  thickness, 
isotropic  model  with  clamped  boundary  conditions  at  the  lower 
boundary.  The  model  uses  primarily  QUAF  410  plate  elements  to 
approximate  the  surface  with  TRINC  320  triangular  elements  (see 
reference  [16])  used  where  necessary. 

The  constant  thickness  model  actually  has  two  different 
thicknesses  associated  with  it.  The  row  of  elements  along  the  entire 
lower  edge  are  0.29  in.  thick.  The  elements  in  the  rest  of  the  shell 
are  0.1455  in.  thick. 

The  variable  thickness  model  better  approximates  the  actual 
thickness  distribution  in  the  shell.  Figure  3.12  shows  this 
thickness  variation  over  the  shell.  A  taper  was  built  into  the 
actual  composite  shell  between  thickness  regions  to  gradually  change 


50 


51 


Figure  3.12  Shell  Thickness  Variation 


the  thickness  variations;  eliminating  stress  concentrations  due  to 
step  transitions.  To  approximate  this  thickness  variation,  without 
refining  the  grid,  an  average  element  thickness  was  used  based  upon 
the  element  surface  area.  On  a  scale  drawing  similar  to  Figure  3.12, 
elements  that  had  a  thickness  variation  line  cutting  through  them 
were  given  average  thicknesses .  For  example ,  if  an  element  was  cut 
so  that  60%  of  the  element  area  was  in  the  0.1455  in.  area  and  the 
remaining  40%  was  in  the  0.1419  in  area  the  resulting  thickness  would 
be : 

t  -  0.6(0.1455  in.)  +  0.4(0.1419  in.)  -  0.144  in. 

avg 

Thicknesses  were  rounded  to  three  decimal  points  except  in  areas 
where  they  are  completely  contained  in  a  thickness  region.  In  this 
case  they  were  given  the  associated  four  decimal  point  thickness 
value.  The  thickness  variation  modeling  resulted  in  26  different 
thicknesses;  11  at  the  lower  edge  and  15  in  the  upper  section  of  the 
model . 

As  stated  earlier,  the  models  in  this  phase  were  developed  to 
better  understand  the  STAGS  program  and  the  analysis  in  general. 
Therefore,  all  models  in  Phase  I  use  an  isotropic  material  modeling 
(orthotropic  modeling  is  considered  in  sections  3.3  and  3.4)  to 
simplify  the  analysis. 

The  material  properties  for  the  Kevlar/Polyester  composite 
fabric  were  determined  through  testing  by  the  Air  Force  Wright 
Aeronautics  Laboratory  (AFWAL)  Materials  Laboratory,  Wright-Patterson 
AFB,  Ohio.  The  results  of  tensile  tests  on  the  material  gave: 

E  -  4.89x10®  psi 

G  -  0.19x10®  psi 
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1/  -  0.038 


Since  the  Kevlar/Polyester  composite  weave  is  not  an  isotropic 
material,  some  assumptions  were  made.  The  value  for  Young's  modulus, 
E,  was  used  as  shown.  Poisson's  ratio,  i/,  was  modified  to 
approximate  an  isotropic  material  such  as  aluminum  with  a  value  of 
0.33  [19].  The  shear  modulus,  G,  was  also  modified  to  fit  its 
isotropic  definition  of  [14] 


2(1+^/) 


(3.1) 


Using  the  value  for  E  and  the  assumed  u  value  gives,  from  Equation 
(3.1),  a  shear  modulus  of  2 . 362xl0®psi . 

As  shown  in  Table  3.2,  the  boundary  conditions  used  in  this 
phase  of  modeling  was  a  clamped  boundary  condition  around  the  entire 
bottom  edge.  Clamped  means  that  all  boundary  degrees  of  freedom  are 
constrained  to  zero  displacement. 

The  QUAF  410  quadrilateral  element  and  the  TRINC  320  triangular 
element  were  used  in  this  Phase  due  to  their  reduced  degrees  of 
freedom  as  compared  to  the  QUAF  411  and  TRINC  321  elements.  Also, 
with  the  new  rotational  degree  of  freedom  formulation  (see  section 
2.2.5)  the  410  element  is  better  behaved  under  rotation  [18]. 

The  Phase  I  models  were  run  on  the  VAX  11/785  computer  using  the 
STAGSC-1  program.  One  immediate  result  from  the  runs  was  the  amount 
of  computer  time  (Central  Processing  Unit  (CPU)  time)  necessary  to 
complete  these  runs.  The  linear  run  took  only  11  minutes  as  opposed 
to  the  nonlinear  runs  which  took  4  hours  and  37  minutes  each.  As 
stated  earlier,  the  nonlinear  runs  used  the  full  Newton- Raphson 
method  (as  opposed  to  the  modified  Newton- Raphson  method)  as  outlined 
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in  section  2.2.5. 


Load  versus  normal  displacement  plots  are  shown  in  Figures  3 . 13 
and  3.14  for  nodes  168  and  176  (Figure  3.5)  respectively. 

From  examining  Figures  3.13  and  3.14  it  is  evident  that  there  is 
a  difference  between  the  linear  and  nonlinear  constant  thickness 
models.  The  nonlinear  model  is  stiffer  due  to  the  coupling  of 
in-plane  and  bending  forces  due  to  the  higher  order  terms  in  the 
nonlinear  kinematic  relations. 

The  nonlinear  models  also  show  a  difference  due  to  the  thickness 
variation.  At  the  particular  nodes  graphed  in  Figures  3.13  and  3.14, 
the  thickness  difference  is  0.0144  in.  This  effect  is  not  as 
dramatic  as  the  previously  mentioned  effect  of  a  linear  versus 
nonlinear  run,  but  it  does  show  a  reduction  in  the  stiffness  of  the 
model  in  this  area  due  to  the  reduction  of  the  bending  moments ; 
caused  by  the  thinner  model.  There  also  is  an  interesting  difference 
between  the  linear  model  and  the  nonlinear  variable  thickness  model 
below  about  0.4  of  the  total  load.  Initially  the  nonlinear  model  is 
not  as  stiff  as  the  linear  model,  but  this  changes  after  about  0.4  of 
the  total  load  due  to  the  nonlinear  coupling  of  membrane  and  bending 
forces . 
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LACEMENT  (IN.) 


FRACTION  OF  TOTAL  LOAD 


Figure  3.14  Node  176  Phase  I  Modeling;  Isotropic  Material, 
410  Element  Clamped  B.C.s 


3 . 3  Phase  II  Modeling 


The  basic  model  carried  forward  from  the  Phase  I  modeling  is  the 
nonlinear  variable  thickness  model.  All  of  the  models  in  Phase  I 
were  modeled  with  an  isotropic  assumption,  clamped  boundary 
conditions,  and  QUAF  410  quadrilateral  elements.  The  basic  changes 
to  this  model  include  modifying  the  isotropic  material  assumption 
(modeling  it  as  an  orthotropic  material)  and  modifying  the  boundary 
conditions.  The  three  models  developed  in  Phase  II  are  shown  in 
Table  3.3 


Mod 

Solution 

Thickness 

Material 

B.C.s 

Elem 

4 

Nonlin. 

Variable 

Orthotrop 

Clamp 

410 

2888(1908) 

362 

5 

Linear 

Variable 

Orthotrop 

Clamp 

410 

2888(1908) 

362 

6 

Nonlin. 

Variable 

Orthotrop 

Mod. 

410 

3616(2268) 

466 

*D0F«  Degrees  of  Freedom  Elem*  Total  #  of  Elements 

ADOF*  Active  Degrees  of  Freedom  Type*  Element  Type 


Table  3.3  Phase  II  Models 


The  actual  material  used  in  the  shell  construction  is  a 
Kevlar-49/F- 141  Polyester  composite  fabric  (weave).  The  edge  (0.29m 
thickness  area  in  Figure  3.12)  also  includes  a  one  ply  inner  and 
outer  edge  backup  of  Glass  7781/F-141  Polyester.  STAGS  is  capable  of 
handling  composite  materials  through  the  thickness  of  the  shell  (the 
wall  composition) ,  but  the  specific  types  include  a  layered  wall 
(orthotropic),  a  fiber  wound  wall,  and  a  corrugated  wall  [16],  The 
shell  material  is  a  layered  orthotropic  material  (composite  weave) 
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but  this  Is  not  the  type  of  orthotropic  material  that  STAGS  is 
designed  to  handle  (STAGS  is  designed  for  unidirectional  composite 
materials) . 

After  a  thorough  search  on  how  to  model  this  weave  material,  two 
solutions  were  found.  The  first  was  to  model  the  weave  as  a  layered 
isotropic  material  using  all  three  of  the  material  constants  (E,  G, 
and  i/)  calculated  by  the  AFWAL  Materials  Laboratory  (see  section 
3.2).  The  second  was  to  model  the  structure  as  a  unidirectional 
orthotropic  layered  material  which  would  involve  modifying  the  ply 
lay-up  and  estimating  a  few  material  constants. 

After  examining  a  few  test  cases  of  modeling  a  simple  structure 
with  the  same  types  of  weave  solutions  as  previously  mentioned,  it 
was  found  that  the  orthotropic  (unidirectional  and  layered)  material 
modeling  was  much  more  flexible  than  the  isotropic  modeling. 

With  the  test  cases  in  mind  and  from  looking  at  preliminary 
experimental  data  in  the  same  region  (showing  very  large  deflections) 
the  orthotropic  material  modeling  was  used  for  the  Kevlar/Polyester 
material . 

To  model  the  Kevlar/Polyester  composite  orthotropically ,  each 
lamina  (ply)  was  split  in  half  and  modeled  as  two  0  and  90  degree 
plies  (each  having  one-half  the  thickness  of  the  original  weave  ply) . 
At  first  glance  this  may  not  seem  like  a  good  approach  since  the 
lay-up  will  not  be  symmetric  (the  inner  ply  will  always  be  at  0 
degrees  orientation  and  the  outer  ply  at  90  degrees)  and  a  coupling 
of  in-plane  and  bending  forces  due  to  the  material  matrix  [B]  (see 
Appendix  C)  will  result.  This  effect  occurs  but  is  very  mild  due  to 
the  number  of  plies  in  the  lay-up  and  the  thickness  of  the  plies.  In 


58 


order  to  model  the  shell  this  way  the  properties  for  an  orthotropic 
material  were  found.  For  an  orthotropic  material  (plane  stress 
condition)  the  properties  are  [14]: 

-  Young's  modulus  in  the  1-direction 

-  Young's  modulus  in  the  2-direction 

V  -  Poisson's  ratio  in  the  2-direction  when  loaded  in 
12 

the  1- direction 

u  -  Poisson's  ratio  in  the  1-direction  when  loaded  in 
21 

the  2-direction 

G  -  Shear  modulus  in  the  1-2  plane 

where  the  subscripts  1  and  2  indicate  material  coordinates  as  shown 
in  Figure  C.2,  Appendix  C.  In  order  to  find  these  values,  some 
approximations  were  made. 

From  the  Phase  I  modeling  section,  the  results  of  the  AFWAL 
tests  on  the  Kevlar/Polyester  composite  were: 

E  —  4.89x10^  psi 
G  -  0.19x10®  psi 
1/  -  0.038 

The  Young's  modulus,  E,  from  above  can  be  assumed  as  E^  since  it  is 

the  modulus  in  the  fiber  direction  and  for  a  weave  this  is  0  or  90 

degrees.  From  examining  Poisson's  ratios  for  composite  materials  in 

references  [14]  and  [20]  it  would  appear  that  the  Poisson's  ratio 

given  would  correspond  to  due  to  its  extremely  low  value. 

Therefore  the  given  Poisson's  ratio  will  be  used  as  The  shear 

modulus,  G,  is  the  in-plane  shear  therefore  it  will  be  used  as  G  . 

The  only  value  left  to  find  is  E^;  since  for  an  orthotropic  material 

1/  can  be  found  from  [14] 

12  ^  ' 
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(3.2) 


Orthotropic  material  properties  could  not  be  found  for 
Kevlar /Polyester ,  however  Graphite/Epoxy  values  are  readily 
available.  To  estimate  a  value  of  for  Kevlar/Polyester  a  ratio  of 
E^  to  Graphite/Epoxy  was  used  along  with  the  assumed 

Kevlar/Polyester  E^  value.  From  reference  [20]  a  ratio  of  E^  to  E^ 
for  Graphite/Epoxy  is  approximately 


With  the  ratio  found  in  Equation  (3.3)  and  the  Kevlar/Polyester 

value,  E  becomes 
2 


E 

2 


E 

_i _ 

11.56 


4.89x10 

11.56 


4.23x10^ 


(3.4) 


Now  if  one  uses  Equation  (3.2)  a  value  for  can  be  found  resulting 
in 


Properties  for  the  edge  backup  material  (Glass/Polyester  cloth) 
used  on  the  shell  were  not  determined  experimentally,  therefore 
approximate  values  were  used. 

The  glass  plies  were  treated  as  an  isotropic  material,  as  in 
Phase  I  for  Kevlar/Polyester,  since  they  only  make  up  about  six 
percent  of  the  edge  thickness.  A  value  for  Young's  modulus  for 
Glass/Polyester  cloth  was  found  in  reference  [21]  and  was  3.24x10^ 
psi.  A  Poisson's  ratio  could  not  be  found  for  a  cloth  material  so  a 
Poisson's  ratio  was  estimated,  as  in  Phase  I  modeling,  as  0.33.  The 
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shear  modulus  was  calculated  from  Equation  (3.1)  as  1.22x10^  psi. 

In  order  to  model  the  shell  as  a  laminate  of  composite  plies  the 
specification  for  the  ply  lay-up  was  used.  A  listing  for  the  ply 
lay-up  is  shown  in  Table  3.4  for  the  different  thickness  areas  shown 
in  Figure  3.12.  This  table  is  for  the  actual  shell;  the  model  is 


SIEflR 

total 

t 

Mat. 

Region  Thicknesses  (Figure  3.12)  J 

0.1383 

0.1455 

0.2900 

At 

.0085 

.0085 

G 

— 

X 

1 

.0036 

.0036 

K 

X 

X 

X 

X 

X 

2-16 

.0085 

.1275 

K 

X 

X 

X 

X 

X 

17-18 

.0036 

.0072 

K 

X 

X 

X 

X 

19 

.0036 

.0036 

K 

X 

X 

X 

20-34i 

.0085 

.1275 

K 

X 

35 

.0036 

.0036 

K 

X 

X 

_ B± _ 

.0085 

.0085 

G 

X 

Total 

Plies 

16 

18 

19 

20 

37 

♦ 

K-Kevlar/Polyester  f  Numbered  from  outer  ply  inward 
G«Glass/Polyester  J  Used  only  on  the  edge 
(Dimensions  in  inches) 


Table  3.4  Actual  Shell  Ply  Lay  Up 


layed  up  the  same  way  except  that  each  ply  in  the  finite  element 
model  is  made  up  of  two  0  and  90  degree  plies  each  having  half  the 
thickness  of  the  actual  shell  ply. 

As  an  example  of  the  lay-up  used  in  the  finite  element  model 
consider  an  element  in  the  0.1419  in.  thickness  region  (see  Figure 
3.12).  Doubling  the  number  of  plies  and  halving  each  ply  thickness 
in  Table  3.4  gives  the  lay-up  in  this  thickness  region,  for  the 
finite  element  model,  as  follows: 

-  2  layers  of  0.0018  in.  Kevlar/Polyester 
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-  30  layers  of  0.00425  in.  Kevlar/Polyester 

-  4  layers  of  0.0018  in.  Kevlar/Polyester 

-  2  layers  of  0.0018  in.  Kevlar/Polyester 

where  all  plies  alternate  between  0  and  90  degrees  orientation. 

Thicknesses  that  lie  between  the  actual  shell  thicknesses  (see 
Figure  3.12)  due  to  the  variable  thickness  approximation  (see  section 
3.2)  are  modeled  as  in  the  following  example  of  a  hypothetical  0.144 
in.  thick  element.  This  thickness  lies  between  the  0.1419  in.  and 
0.1455  in.  thickness  regions  in  Table  3.4.  The  model  is  first  given 
the  total  lay-up  of  the  0.1419  in.  region  as  in  the  previous  example. 
The  remaining  thickness  (0.0021  in.)  is  divided  into  two  0  and  90 
degree  plies  (corresponding  to  ply  no.  35,  Table  3.4)  of 
Kevlar/Polyester . 

All  thicknesses  above  0.1455  in.  in  the  variable  thickness  model 
correspond  to  edge  thicknesses.  These  approximated  thicknesses  for 
the  edge  elements  vary  from  0.231  to  0.29  in.  From  Table  3.4  the 
total  number  of  plies  in  the  actual  shell  edge  is  37  which 
corresponds  to  74  plies  in  the  orthotropic  approximation  for  the 
model.  After  trying  a  74  ply  run  on  the  computer  it  was  found  after 
several  attempts  that  STAGS  is  limited  to  approximately  50  plies. 

With  this  in  mind  the  edges  were  modeled  with  42  plies.  Plies  1, 
17-18,  and  35  were  left  as  shown  in  Table  3.4  (modified  as  before  by 
halving  each  shell  ply) .  Plies  A  and  B  were  left  at  full  thickness 
and  treated  isotropically.  Plies  2-16  in  Table  3.4  were  modeled  as 
20  plies  at  0.006375  in.  each.  Plies  20-34  in  Table  3.4  are  where 
the  variable  thicknesses  on  the  edge  were  taken  into  account.  They 
were  modeled  as  10  plies  with  variable  thicknesses  to  accommodate 


62 


thicknesses  on  the  edge  other  than  0.29  in. 

All  of  the  models  in  Phase  I  and  models  4  and  5  in  Phase  II  have 
clamped  boundary  conditions  at  the  lower  edge.  A  cross  section  of 
the  actual  boundary  (continuous  around  the  entire  lower  edge)  on  the 
lower  edge  of  the  shell  is  shown  in  Figure  3.15.  The  lower  edge  of 
the  finite  element  model  is  shown  in  the  figure  and  this  is  where  the 
clamped  models  (Models  1-5)  were  assumed  fixed.  The  angle  d  shown  in 
Figure  3.15  is  zero  degrees  over  much  of  the  boundary  except  at  each 
end  of  the  shell  where  it  approaches  ten  degrees. 

By  examining  Figure  3.15  it  is  obvious  that  the  clamped  boundary 
condition  that  was  originally  assumed  is  too  stiff.  The  0.5  in. 
aluminum  plate  will  allow  lateral  deflection  due  to  bending,  and  the 
rubber  washer  in  the  plane  of  the  shell  will  allow  vertical 
deflection.  These  two  components  are  changed  in  the  modified 


Figure  3.15  Shell  Lower  Boundary 


boundary  conditions. 

To  model  the  aluminum  plate,  a  0.5  in.  high  by  0.125  in.  thick 
plate  element  was  added  to  the  finite  element  model  at  the 
lower  edge  (see  Figure  3.16).  The  lengths  of  the  elements  were 
determined  by  the  original  lower  boundary  elements  directly  above  the 
aluminum  plate  elements.  Degrees  of  freedom  at  the  interface  of  the 
shell  and  plate  were  assumed  free.  The  material  properties  of  the 
aluminum  plate  are  [22] 

E  -  10. 7x10®  psi 

1/  -  0.33 

G  -  4.0x10®  psi 

To  simulate  the  rubber  washer  a  STAGS  general  beam  element 
(GSBM2  220)  was  used.  These  beam  elements  were  attached  to  the  lower 
nodes  of  the  aluminum  plate  element  and  are  oriented  vertically  (see 
Figure  3.16).  Since  material  properties  for  the  actual  rubber  could 
not  be  found  and  even  if  they  could  be  found,  an  analytic  analysis  of 
the  equivalent  stiffness  of  the  washers  would  be  extremely  difficult. 
Therefore,  an  estimate  had  to  be  made  for  the  equivalent  beam 
stiffness  necessary  to  model  the  washers. 

This  approximation  was  done  using  the  force  (F),  stiffness  (k) , 
and  displacement  (x)  relation  for  a  spring  as  in 

F  -  kx  (3.6) 

If  the  displacements  due  to  the  rubber  washers  could  be  found  in  the 
actual  model  due  to  the  applied  forces,  the  stiffness  could  be  found. 
This  stiffness  would  then  have  to  be  related  to  the  equivalent 
stiffness  in  a  linear  beam  as  in 
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(3.7) 


where  A  is  the  cross-sectional  area,  E  is  the  Young's  modulus,  and  L 
is  the  length  of  the  beam. 

During  experimental  testing  of  the  shell  (see  section  4)  rough 
measurements  were  made  at  the  lower  edge  of  the  shell  to  see  how  much 
vertical  displacement  was  caused  by  the  rubber  washers .  Vertical 
equilibrium  forces  were  then  found  from  the  linear  model  5  (see  Table 
3.3),  using  the  STAGS  post  processor,  at  the  locations  where  the 
rough  measurements  were  taken  experimentally.  Using  these 
experimental  displacements  and  the  finite  element  forces  an 
approximation  of  the  stiffnesses  in  the  rubber  washers  could  be  found 
from  Equation  (3.6).  The  average  stiffness  found  from  Equation  (3.6) 
was  6,292  lb. /in. 

The  beam  used  in  the  model  could  then  be  designed  with  any 
combination  of  variables  shown  on  the  right  hand  side  of  Equation 
(3.7)  equaling  the  estimated  rubber  washer  stiffness.  The  values 
used  for  the  "rubber  beam"  in  the  analysis  were 

A  -  1.15  in^ 

E  ”  2.07x10*"  psi 

L  -  0.377  in. 

Finally,  the  beam  was  modeled  in  STAGS  with  a  linear  kinematic 
relation  since  the  stiffnesses  were  calculated  linearly.  The  beam  is 
clamped  at  the  bottom  and  only  allowed  vertical  deflection  at  the 
top;  to  simulate  the  vertical  deflection  of  the  washer. 

Figure  3.16  shows  the  modified  finite  element  boundary  condition 
in  a  side  and  cross  section  view.  The  figure  shows  the  relation  in 
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the  model  between  the  shell,  the  aluminum  plate,  and  the  "rubber 
beam."  The  additional  elements  (aluminum  plate  and  "rubber  beam")  do 
not  have  an  applied  load;  the  shell  is  loaded  as  in  the  clamped 
models . 


Side  View 


Aluminum  Plate 


0, 


f 


She 


Original  Boun  dar y 
-Orig.  Model  Clamped 
1  -Mod.  Model  Free 


0.5  in 


/ 


Vertical  Deflection  Only 


Rubber  Beam  Boundary 


i  e  w 


Figure  3.16  Modified  Finite  Element  Boundary  Condition 


66 


As  in  Phase  I  modeling,  all  models  in  this  Phase  were  ran  on  the 
VAX  11/785  computer.  The  CPU  time  for  the  nonlinear  clamped  model 
(Model  4)  was  about  5  hours,  29  minutes.  The  same  model  in  a  linear 
run  (Model  5)  took  about  15  minutes.  The  orthotropic  model  with  the 
modified  boundary  conditions  (Model  6)  took  about  6  hours,  34 
minutes . 

Load  versus  normal  displacement  plots  are  shown  in  Figures  3.17 
and  3.18  for  nodes  168  and  176  (Figure  3.5)  respectively. 

From  examining  Figures  3.17  and  3.18  there  is  more  than  a  100 
percent  difference  in  the  linear  versus  nonlinear  runs  for  the 
clamped  models.  The  linear  run  is  assumed  in  error  due  to  its 
displacement  as  compared  to  the  thickness  of  the  shell  (0.1311  in.) 
in  the  area  surrounding  the  node  points  168  and  176.  Linear  plate 
(STAGS  "shell"  element)  theory  assumes  that  the  displacements  are 
much  smaller  than  the  thickness  of  the  shell  [1].  For  this 
particular  area  the  opposite  is  true;  the  displacements  are  much 
larger  than  the  shell  thickness.  Because  of  this  violation  in  the 
thin  plate  assumption  for  linear  theory  ,  and  the  tendency  in  this 
shell  for  large  displacements  under  the  given  loading,  Model  5  was 
the  last  linear  run  attempted. 

The  difference  between  the  clamped  and  modified  boundary 
condition  runs  (Figure  3.17  and  3.18)  shows  a  slight  but  pronounced 
reduction  of  stiffness  in  the  modified  boundary  condition  model.  The 
clamped  boundary  condition  obviously  is  stiffer. 

Finally,  by  going  back  to  Phase  I  and  examining  Figures  3.13  and 
3.14  for  the  same  nodes,  there  is  a  dramatic  difference  between 
modeling  the  structure  isotropically  versus  orthotropically .  The 
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Figure  3.17  Node  168  Phase  II  Modeling:  Orthotropic  Material 
Variable  Thickness,  410  Element 
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Figure  3.18  Node  176  Phase  II  Modeling:  Orthotropic  Material 
Variable  Thickness,  410  Element 


displacements  in  the  orthotropic  model  nearly  double  those  of  the 
isotropic  model.  From  this  comparison  it  is  a  possibility  that  the 
material  assumption  used  in  Phase  I  is  much  too  stiff  since  only  one 
of  the  actual  material  constants  (Young's  Modulus)  was  actually  used. 
In  Phase  II  all  of  the  actual  material  constants  were  used  and  two 
were  assumed  from  the  relation  of  another  orthotropic  material 
(Graphite/Epoxy) .  Comparisons  of  finite  element  solutions  versus 
experimental  results  (section  5)  will  show  if  the  orthotropic 
assumption  is  valid. 


3.4  Phase  III  Modeling 

In  this  phase  the  sixth  model  (Phase  II)  in  the  analysis  is  used 
for  comparison  to  the  two  models  developed  in  this  stage.  Both  of 
the  models  developed  in  Phase  III  are  similar  to  Model  6  with  the 
exception  that  one  model  is  changed  to  model  the  Kevlar/Polyester 
material  as  a  layered  isotropic  material  (to  compare  against  the 
orthotropic  material  model  of  section  3.3)  and  the  other  model 
changes  the  STAGS  element  types.  A  summary  of  these  finite  element 
models  is  shown  in  Table  3.5  along  with  Model  6  from  Phase  II 
modeling. 

The  layered  isotropic  model  (Model  7)  uses  the  same  lay-up  as  in 
Model  6.  The  only  change  to  the  model  is  that  all  of  the  material 
constants  calculated  by  the  AFWAL  Materials  Laboratory  are  used. 

They  are  written  here  again  for  convenience  as 

E  -  4 . 89x10^  psi 
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G  -  0.19x10^  psi 
»/  -  0.038 


Mod 

Solution 

Thickness 

Material 

B.C.s 

ISH 

DOF(ADOF)* 

Elem 

6 

Nonlin. 

Variable 

Ortho trop 

Mod. 

410 

3616(2268) 

466 

7 

Nonlin. 

Variable 

Mod. 

410 

3616(2268) 

466 

8 

Nonlin. 

Ortho trop 

Mod. 

411 

5899(3506) 

466 

*DQF“  Degrees  of  Freedom  £lsm"  Total  #  of  Elements 


ADOF*  Active  Degrees  of  Freedom  Type*  Element  Type 


Table  3.5  Phase  III  Models 


The  ply  orientation  in  the  model  is  left  at  0  and  90  degrees  since  in 
an  isotropic  material  the  material  constants  are  the  same  in  either 
direction. 

The  final  model  in  the  analysis,  Model  8,  is  the  same  as  Model  6 
with  the  exception  that  the  QUAF  410  quadrilateral  element  (see 
section  2.2.4)  is  replaced  by  the  QUAF  411  element,  the  TRINC  320 
triangular  element  with  the  TRINC  321,  and  the  GSBM2  220  beam  element 
with  the  GSBM2  221. 

Models  7  and  8  were  ran  on  the  VAX  11/785  computer  in  the  same 
manner  as  the  other  runs.  The  CPU  time  for  the  layered  isotropic 
model  (Model  7)  was  6  hours  and  9  minutes  which  was  close  to  the  6 
hours  and  43  minutes  it  took  to  run  Model  6.  A  dramatic  difference 
in  CPU  time  was  found  in  Model  8.  Model  8  took  26  hours  and  51 
minutes  to  run.  It  was  anticipated  that  the  Model  8  run  would  take 
longer  due  to  the  increased  number  of  degrees  of  freedom  (see  Table 
3,5),  but  not  as  long  as  it  actually  did. 
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Load  versus  normal  displacement  plots  are  shown  in  Figures  3.19 
and  3.20  for  nodes  168  and  176  (see  Figure  3.5)  respectively. 

From  examining  Figures  3 . 19  and  3.20  it  is  evident  that  the 
layered  isotropic  material  model  (Model  7)  is  a  much  stiffer  model 
than  the  orthotropic  model  (Model  6) .  From  examining  experimental 
data  it  is  also  evident  that  the  orthotropic  model  gives  better 
results  in  terms  of  displacements. 

Finally  from  Figures  3.19  and  3.20  the  change  to  a  higher  degree 
of  freedom  element  (QUAF  411  Element)  was  not  necessary.  In  fact, 
given  the  CPU  time  necessary  to  complete  the  Model  8  (411  Element) 
run  versus  the  time  for  the  Model  6  (410  Element) ,  the  QUAF  411 
element  seems  to  be  very  impractical  for  mild  nonlinearities  (as 
opposed  to  a  shell  collapse  which  is  highly  nonlinear) . 

Since  the  QUAF  410  Element  model  (Model  6)  was  much  more 
economical  and  gave  essentially  the  same  results  as  the  QUAF  411 
model  (Model  8)  and  since  the  displacement  results  for  Model  6  are 
much  better  than  those  for  Model  7,  Model  6  was  used  for  comparison 
with  the  experimental  results  (see  section  5).  The  STAGS  input  deck 
for  Model  6  is  shown  in  Appendix  G.  The  UPRESS  loading  subroutine 
used  to  apply  the  pressure  loading  in  this  model  is  shown  in 
Appendix  H. 
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Figure  3.19  Node  168  Phase  III  Modeling:  Variable  Thickness 
Modified  B.C.s 
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Figure  3.20  Node  176  Phase  III  Modeling:  Variable  Thickness, 
Modified  B.C.s 


The  experimentation  on  the  actual  shell  was  done  by  theAir  Force 
Wright  Aeronautics  Laboratory,  Wright- Patterson  AFB,  Ohio.  The  test 
simulated  the  actual  aerodynamic  loading  on  the  shell  by  using 
inflatable  bladders  within  the  shell  to  provide  internal  pressure 
loading. 

The  steps  necessary  to  devise  a  test  fixture  that  can  produce 
this  type  of  pressure  loading  is  as  follows.  First,  an  aluminum 
frame  (Figure  4.1)  was  built  that  was  the  same  shape  as  the  shell  but 


was  smaller  so  that  the  shell  would  fit  over  the  frame  and  leave 


about  1  1/2  in.  between  the  two.  The  frame  also  had  spaces  in  it 
that  coincide  with  the  load  regions  shown  in  Figure  3.7.  Next,  holes 
were  drilled  into  each  load  area  of  the  test  frame.  Rubber  hoses, 
extending  from  electronically  controlled  pressure  valves,  were 
inserted  into  each  hole  and  the  gap  around  the  hoses  sealed.  This 
entire  assembly  was  then  bolted  to  a  test  stand. 

Rubber  sheets  were  then  glued  to  the  frame  with  contact  cement; 
covering  each  load  region  and  forming  the  basic  loading  bladder 
(Figure  4.2).  A  border  of  one  inch  thick  closed  cell  foam  was 
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Figure  4.2  Rubber  Bladders 


attached  (glued  to  the  rubber  and  tied  to  the  aluminum  frame)  to 
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the  intersection  of  each  load  region.  This  foam  border  was  about 
twelve  inches  wide  (extending  six  inches  into  each  region)  tapered 
from  one  inch  thick  at  its  midpoint  (the  load  region  boundary)  to  a 
zero  thickness  at  the  edges.  This  foam  border  kept  the  bladder  in 
one  region  from  bulging  into  another  region.  The  taper  kept  the 
inner  edges  of  the  closed  cell  foam  from  touching  the  shell  before 
the  bladder  contacts  the  shell. 

The  shell  was  then  placed  over  the  frame  and  bolted  to  the  test 
stand  (Figure  4.3)  with  same  type  of  boundary  restraint  as  shown  in 


Figure  4.3  Experimental  Set-Up 


Figure  3 . 15.  Finally ,  linear  variable  differential  transducers  (LVDTs) 
were  placed  around  the  mid-section  of  the  shell  (Figure  4.4).  The 
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Figure  4.4  Close-Up  of  Displacement  Transducers 

shell  at  specified  points  that  coincide  with  node  points  on  the 
finite  element  model. 

During  loading  it  was  anticipated  that  the  bladders  would  tend 
to  bulge  more  toward  the  center  of  a  load  region;  reducing  the 
contact  area  of  the  pressure  load.  To  compensate  for  this,  the  4950 
Test  Wing  recalculated  the  loading  (Table  3.1)  on  the  shell  to 
account  for  the  reduced  contact  area  of  the  bladder  as  compared  to 
the  true  pressure  values.  Table  4.1  shows  modified  loads  used  in  the 
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finite  element  modeling  (Figure  3.5). 
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Table  4.1  Modified  Static  Equivalent  Loads 

A  computer  was  used  to  control  the  loading  on  the  shell .  The 
computer  would  open  and  close  individual  valves  that  controlled  the 
air  pressure  to  each  load  region.  Pressure  transducers  provided 
feedback  to  the  computer  to  insure  that  the  regions  were  loaded 
properly. 

During  testing  the  shell  was  loaded  in  increments  (10%  of  total 
load)  up  to  the  total  load  and  back  down  to  an  unloaded  condition. 

The  computer  that  controlled  the  load  on  the  shell  also  took  readings 
from  all  of  the  measuring  devices  on  the  shell  at  predetermined  time 
steps.  Displacement  results  obtained  during  testing  are  shown  along 
with  finite  element  solutions  in  the  following  section. 
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As  s  ted  in  the  Phase  III  Modeling  section  (section  3.4),  the 
finite  element  model  used  for  comparison  to  experimental  results  is 


the  variable  thickness,  orthotropic  material  model  with  modified 
boundary  conditions  (Model  6)  and  a  QUAF  410  element.  In  order  to 
show  a  general  deformation  pattern  for  the  shell  a  cross-section  view 
of  the  undeformed  and  deformed  shape  of  the  shell  is  shown  at  the 
same  cross-section  (constant  x  value)  as  the  nodes  that  were  shown  in 
the  Phase  I-III  modeling  sections  (see  Figure  3.5). 

From  Figure  5.1  a  general  deformation  pattern  is  apparent.  The 
shell  is  bulging  at  the  sides  and  the  top  is  deflecting  down.  Also 
from  examining  the  data  from  the  model  and  experimentation,  the 
displacements  on  the  sides  of  the  shell  at  a  y  value  of  approximately 
eleven  inches  are  the  largest.  These  displacements  are  an  order  of 
magnitude  larger  than  the  displacements  at  the  top  of  the  shell. 

There  also  seems  to  be  a  lack  of  consistency  in  the  experimental 
deflections  at  the  top  of  the  shell.  Because  of  the  large  difference 
in  displacements  and  the  inconsistency  in  the  experimental 
displacements  at  the  top  of  the  shell,  the  side  node  displacements 
will  be  compared.  Finally,  a  node  on  the  top  of  the  shell  is  shown 
for  general  trends  in  this  area. 

Comparison  of  load  versus  displacement  are  made  for  the 
nodes  shown  in  Figure  5.2.  The  load  versus  displacement  curves  for 
these  nodes  are  shown  in  Figures  5.3  to  5.9. 
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Figure  5.1  Deformed  vs.  Undeformed  Finite  Element  Cross  Section 

From  examination  of  Figures  5.3  to  5.8  for  the  side  nodes,  the 
finite  element  solutions  are  within  12.5%  of  the  experimental  results 
at  maximum  load.  The  Figures  also  show  a  "hysteresis"  type  effect 
for  the  experimental  results.  This  effect  is  not  a  true  hysteresis 
(plastic  deformation) ;  it  is  caused  by  the  shell  moving  at  the 
boundaries  because  of  the  rubber  washers  (see  section  3.4)  and  not 
settling  completely  back  into  the  boundary  attachment. 

These  nodes  also  show  very  large  displacements  (greater  than  0.5 


8 


in.)  which  are  much  larger  than  the  thickness  of  the  shell.  For 


163  219(227) 


Figure  5.2  Shell  and  Finite  Element  Nodes  Used  for  Comparison 

example  the  displacement  at  node  168  is  greater  than  0.6  in.  and  the 
thickness  is  only  0.1311  in.  This  type  of  displacement  justifies  the 
use  of  a  nonlinear  finite  element  model  since  the  basic  assumption 
behind  a  linear  analysis  is  that  the  displacements  must  be  much 
smaller  than  the  thickness  of  the  shell  [Ij. 
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Figure  5.3  Node  151:  Experimental  vs.  Finite  Element 
Surface  Normal  Displacement 
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Figure  5.4  Node  159:  Experimental  /s .  Finite  Element 
Surface  Normal  Displacement 
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Figure  5.5  Node  168;  Experimental  vs.  Finite  Element 
Surface  Normal  Displacement 
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Figure  5.7  Node  219:  Experimental  vs.  Finite  Element 
Surface  Normal  Displacement 


Figure  5.8  Node  227:  Experimental  vs.  Finite  Element 
Surface  Normal  Displacement 
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From  looking  at  Figure  5.9  for  node  163  on  the  top  of  the  shell 
the  experimental  value  at  maximum  load  is  much  smaller  than  the 
finite  element  solution.  But,  the  experimental  displacements  show  a 
very  erratic  behavior  at  this  and  other  top  node  locations .  The  only 
part  of  this  loading  curve  that  makes  sense  is  the  section  from 
approximately  0.0  to  0,2  of  the  total  load.  The  curve  shows  an 
upward  movement  of  the  shell  due  to  the  rubber  washer  until  full 
compression  of  the  washer  at  about  0.2  of  the  total  load.  This 
phenomenon  was  noted  during  testing;  a  general  upward  movement  of  the 
shell  at  the  beginning  and  then  movement  of  the  shell  due  to 
kinematics.  This  general  pattern  is  also  evident  in  the  side  nodes 
(Figures  5.3  to  5.8).  Up  until  about  0.1  of  the  total  load  there  is 
no  lateral  displacement  (the  side  node  normal)  at  the  side  nodes. 
Above  this  point  the  side  nodes  start  deflecting  outward  causing  the 
bulging  of  the  shell  as  depicted  in  Figure  5.1.  What  this  also  says 
is  that  the  modeling  of  the  rubber  washers  using  "rubber  beams"  (see 
section  3.4),  that  deflect  during  the  entire  load  process,  will  not 
accurately  model  this  type  of  movement;  essentially  rigid  body 
motion.  Above  about  0.2  of  the  total  load  the  curve  does  not  follow 
a  definite  load  path  and  is  considerably  different  from  that  of  the 
unloading  path.  After  discussion  with  the  Wright  Aeronautics 
Laboratory,  there  does  not  seem  to  be  any  definite  reason  for  this 
load  behavior.  However,  it  must  be  emphasized,  once  again,  that 
these  displacements  are  extremely  small  (in  the  hundredths  of  inches) 
and  are  of  an  order  of  magnitude  less  than  the  side  displacements. 
Therefore,  the  writer  feels  that  they  are  not  important  relative  to 
the  displacements  on  the  side  of  the  shell. 
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NORMAL  DISPLACEMENT  (IN.) 


0.010 


Figure  5.9  Node  163:  Experimental  vs.  Finite  Element 
Surface  Normal  Displacement 
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Overall  the  results  from  the  finite  element  model  using  STAGSC-1 
are  very  good.  They  also  show  that  for  a  shell  as  large  and  flexible 
as  the  one  considered  in  this  thesis  that  a  nonlinear  analysis  is 
essential  in  accurate  modeling. 
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An  important  factor  in  the  post  processing  of  a  f initeelement 
analysis  is  the  determination  of  whether  the  loads  on  the  structure 
are  large  enough  to  cause  failure.  This  is  especially  true  in 
composite  materials  since  they  are  usually  used  in  areas  that  require 
a  tailoring  of  material  properties  in  a  particular  direction  to  meet 
a  load  requirement.  Simply  looking  at  the  stresses  in  the  structure 
from  a  computer  printout  and  comparing  them  to  uniaxial  material 
properties  is  inadequate  due  to,  in  general,  a  biaxial  if  not 
triaxial  stress  state  in  the  actual  structure  [6] .  To  analyze  a 
structure  for  possible  failure  under  these  conditions,  biaxial 
strength  theories  have  been  developed.  One  such  strength  theory  is 
the  Tsai-Wu  theory  that  deals  with  failure  in  an  individual  ply  [6], 
In  the  Tsai-Wu  theory,  failure  in  a  ply  occurs  if  the 
combination  of  strength  tensors  (F^,  etc.)  and  the  stresses  on  the 
ply  in  the  material  coordinates  ,<7^,  and  are  greater  than 

unity.  In  equation  form,  for  an  orthotropic  ply  under  plane  stress 
conditions,  failure  occurs  if  [6] 

Fa  +  Fa  +  F  a^  +  F  a^  +  F  +2F  a  a  >1  <o.l) 

11  22  11  1  22  2  66  12  12  1  2 


where 


t  C 


F  -  -  — t — 
11  X  X 

t  c 


(6.2) 

(6.3) 

(6.4) 
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(6.5) 


F  -  -  — = — 
'22  Y  Y 


1- 


<7  + 


1  1 


X  X 

t  c 


Y  Y 

t  c 


(6.6) 

(6.7) 


and 


X^  -  ultimate  tensile  strength  in  the  1-direction 
X  -  ultimate  compressive  strength  in  the  1-direction 

C 

Y^  -  ultimate  tensile  strength  in  the  2 -direction 
Y  -  ultimate  compressive  strength  in  the  2 -direction 

C 

S  -  ultimate  longitudinal  shear  strength 
To  find  it  is  necessary  to  find  the  biaxial  tensile  failure 
stress,  a.  Since  this  property  is  difficult  to  find  experimentally 
reference  [23 j  suggests  that 


2(X  X  Y  Y  )^'^ 

t  C  t  C 

To  use  this  theory  for  the  Kevlar/Polyester  shell  analyzed,  a 
few  assumptions  need  to  be  made  since  the  shell  was  modeled  as  an 
orthotropic  material  even  though  it  was  a  composite  weave.  The  Air 
Force  Wright  Aeronautics  Laboratory  calculated  ultimate  strengths  for 
the  Kevlar/Polyester  material  in  tension,  compression,  and  shear. 
These  values  were  used  for  X  ,  X  ,  and  S  respectively  and  are  given 

t  C 

as 


X^  -  70,400  psi 

X  -  13,000  psi  (6.9) 

C 

S  -  8,690  psi 

The  other  ultimate  stresses  are  found  using  the  same  assumption  as 
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that  used  in  section  3.3;  the  properties  of  Graphite/Epoxy  are  used 
to  ratio  the  Kevlar/Polyester  fabric  properties . 


From  reference  [20]  ,  the  ratio  of  ultimates  for  Graphite/Epoxy 

are 


130,989  psi 
X  9,290  psi 

C 

197,143  psi 
Y  35,204  psi 

C 


(6.10) 


One  can  now  divide  the  Kevlar/Polyester  ultimate  strengths,  X  and  X 

t  C 

(Equation  6.9),  by  14.1  and  5.6  (Equation  6.10)  respectively,  to  get 

the  estimated  Y  ,  Y  values  fo^*  Kevlar/Polyester.  These  values  are 
t  c 

given  as 


Y^  -  4,993  psi 
X  -  2,321  psi 

e 


(6.11) 


With  all  of  the  Kevlar/Polyester  ultimate  values  defined  (Equations 
(6.9)  and  (6.11)),  all  that  remains  in  the  analysis  is  to  determine 
the  stresses  within  the  laminate  in  the  material  directions. 


Calculating  the  stresses  in  a  ply  in  the  material  direction 
using  STAGS  presented  little  difficulties  since  STAGS  gives  the  user 
the  option  to  print  out  stresses  in  either  the  local  or  material 
directions.  Using  the  same  finite  element  model  as  the  one  used  to 
compare  experimental  versus  finite  element  solutions  (see  section  5), 
the  stresses  were  obtained  in  the  material  directions  for  the  inner 


and  outer  plies;  potentially  the  highest  stresses  due  to  in-plane  and 
bending  effects.  It  was  noted  that  the  stresses  were  given  for  each 
of  the  surfaces  on  the  inner  and  outer  ply.  Therefore,  an  average  of 
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these  values  was  used  as  the  ply  stress  value.  Also,  since  the  inner 
and  outer  ply  on  the  edges  was  a  Glass/Polyester  material  and  was  not 
of  prime  concern  in  this  analysis  and  since  during  testing  it  was 
visually  verified  that  there  was  no  apparent  composite  failure  in  any 
part  of  the  shell,  the  lower  edge  elements  were  omitted  from  the 
analysis , 

A  Fortran  program  was  written  to  enter  each  of  the  average  inner 
and  outer  ply  stresses  and  the  Kevlar/Polyester  ultimate  stress 
values  (Equations  (6.9)  and  (6.11))  into  the  left-hand  side  of 
Equation  (6.1).  A  Tsai-Uu  value  was  then  obtained  for  each  inner  an 
outer  ply  within  an  element.  These  values  were  then  sorted  in  tlie 
program  in  descending  order  along  with  their  respective  ply  and 
element  numbers. 

The  highest  Tsai-Wu  value  obtained  in  the  analysis  was  0.31  and 
is  far  below  the  unity  value  required  for  failure  in  Equation  (6.1). 
The  top  nine  Tsai-Wu  values  (0.22  to  0.31)  are  shown  in  Figure  6.1 
since  they  show  a  definite  pattern. 


Figure  6.1  Tsai-Wu  Composite  Failure  Analysis 
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All  of  these  values  are  on  the  inside  ply  of  the  location  shown  in 
Figure  6.1  This  makes  sense  in  the  analysis  since  the  largest 
deflection  gradient  occurs  on  this  side  of  the  model  as  shown  in  Figure 
5.1.  Although  these  values  are  well  below  unity  (failure  value)  they 
show  a  definite  area  of  concern  in  the  model  if  the  loading  is 
increased. 

Given  the  assumptions  used  in  the  model,  it  would  appear  the 
under  the  loading  in  this  analysis,  failure  of  the  composite  material 
would  not  occur.  But,  the  material  values  are  assumptions,  therefore 
an  investigation  into  the  true  material  values  needs  to  be  done  to 
insure  better  confidence  in  the  failure  analysis. 
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7. 


Based  on  the  analysis  of  a  large,  thin  composite  shell,  with 
asymmetric  loading,  presented  in  this  thesis,  the  following 
conclusions  can  be  made: 

1.  Modeling  of  the  structure  should  include: 

•  Taking  into  account  the  thickness  variations  within  the 
structure.  Modeling  of  the  thickness  variations  in  a 
structure  by  using  an  average  element  thickness  appears  to 
work  well  especially  when  reducing  the  grid  size  to 
accommodate  the  varying  thicknesses  is  impractical. 

•  Accurate  modeling  of  the  material  used  in  the  structure 
whether  orthotropic  or  isotropic.  Modeling  an  orthotropic 
material  as  an  isotropic  (1  layer  through  the  thickness) 
material  does  not  work  well  especially  when  assumptions 
have  to  be  made  about  the  material  properties.  The  layered 
isotropic  material  modeling  was  much  too  stiff  for  this 
analysis.  Modeling  a  weave  composite  by  splitting  the 
individual  plies  in  half  and  orienting  the  plies  at  0  and 
90  degrees  produced  reasonable  results.  Also,  making 
assumptions  about  Kevlar/Polyester  composite  weave  material 
properties  from  Graphite/Epoxy  composite  (orthotropic) 
values  appeared  to  work  well. 

•  Boundary  conditions  on  a  structure  have  a  large  effect 
on  the  analysis.  These  effects  are  especially  true  when 
the  boundary  conditions  are  elastically  dependent.  Proper 
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modeling  of  the  boundary  conditions  on  the  structure 
especially  when  standard  finite  element  boundary  conditions 
(clamped,  pinned,  free,  etc.)  do  not  model  the  actual 
boundary  conditions  accurately. 

2.  Use  of  the  QUAF  410  element  is  more  practical  for  cases 
involving  mild  geometric  nonlinearities.  The  QUAF  411  Element 
took  nearly  four  times  as  long  to  run  (over  26  CPU  hours)  as  the 
QUAF  410  element  with  little  difference  in  results . 

3.  Visual  observation  after  experimentation  indicated  no 
failure  in  the  composite  shell.  This  was  concurred  with 
analytically  with  a  Tsai-Wu  strength  analysis. 

4.  Use  of  a  geometric  nonlinear  analysis  is  essential  to 
problems  involving  displacements  as  large  as  four  times  the 
shell  thickness.  The  linear  analysis  carried  out  over  predicted 
some  deflections  by  more  than  100  percent. 

5.  The  nonlinear  variable  thickness  model  with  modified 
boundary  conditions  and  an  orthotropic  material  assumption 
modeled  the  shell  reasonably  well. 

6.  STAGSC-1  with  its  updated  Lagrangian  formulation  and 
nonlinear  capabilities  does  a  very  good  job  in  analyzing  a 
general  shell  with  large  displacements. 
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Appendix  A:  Derivation  of  the  STAGSC-1  Nonlinear 
Strain  Displacement  Equations 

Since  one  of  the  primary  strengths  of  the  STAGS  computer  program 
is  its  ability  to  analyze  shells  nonlinear ly,  a  derivation  of  the 
nonlinear  kinematics  is  presented. 

To  start  the  derivation  a  few  assumptions  need  to  be  made. 

STAGS  uses  plate  elements  to  model  a  shell  structure.  These  plate 
elements  are  considered  thin  so  that  the  in-plane  displacements,  u 
and  V,  and  the  normal  displacement,  w,  is  a  function  of  only  two 
space  variables  [7],  A  plane  stress  assumption  is  made  whereby  7  , 

XZ 

7  ,  «  ,  and  a  are  considered  zero.  And  finally,  the  Kirchhoff-Love 

hypothesis  applies  for  strains  away  from  the  midplane.  With  these 
assumptions,  the  presentation  can  now  be  started. 

Consider  a  line  segment  as  show  in  Figure  A.l  oriented  in  its 
undeformed  and  deformed  (represented  with  an  asterisk)  state.  From 
Figure  A.l,  it  is  evident  that  the  deformed  state  can  be  related  to 
the  undeformed  state  by  [7] 

* 

X  -  X  u 

(A.l) 

z  -  z  +  w 

and  the  differential  change  in  distance  by  [7] 
dx  -  dx(l  +  u,  ) 

X 

dz*  -  dx  w,^  (A. 2) 

*  2  *  2  ♦  2 

(ds  )  -  (dx  y  +  (dz  ) 

By  substituting  the  first  two  equations  of  Equation  (A. 2)  into  the 
last  equation  it  follows  that  [7] 
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J  -  '  -  .  (w.^) 

Using  the  definition  of  strain  [7] 

* 

-  (ds  -  dx)/dx 

and  then  rearranging  it  to  get  an  alternate  form  [7] 


(A. 3) 


(A. 4) 


£  +  1/2  e  -  1/2 

X  X 


dx 


-  1 


(A. 5) 


Figure  A.l  Line  Element  in  Undeformed  and  Deformed  State  [7] 

Substituting  Equation  (A. 3)  into  Equation  (A. 5)  and  assiiming  the 

2 

strain  is  small  so  that  £  can  be  eliminated  gives  [7] 

£  -  u,  +  l/2(u.^  +  w/)  (A. 6) 

XX  XX 

So  far  the  analysis  for  this  line  element  has  been  only  in  one 
plane.  If  this  line  element  is  in  the  midplane  of  a  plate  element, 


the  rotation  of  this  line  element  about  the  normal  has  to  be 


considered.  If  this  rotation  is  included  in  Equation  (A. 6),  as  in 
Sanders'  nonlinear  shell  equations,  and  rewritten  to  reflect  midplane 
strain,  the  result  is 


€°  -  u,  +  l/2(u,^  +  w.^  +  (A. 7) 

Z  X  z  z 

where  in  Sanders'  equations  [3] 


<t>  -  l/2(u,  -  V,  )  (A. 8) 

y  X 

The  term  that  should  be  included  in  Equation  (A. 7)  to  represent  this 
rotation  in  the  x-direction  about  the  normal,  for  a  flat  plate,  is 
V,  .  This  does  not  imply  the  use  of  Sanders’  equation,  just  the 

X 

importance  of  defining  a  normal  rotation  term  for  flat  plates 
representing  a  shell.  Substituting  this  term  in  Equation  (A. 7) 
results  in  [7]. 


-  u,  +  l/2(u,^  +  w,^  +  v,^) 

XX  XXX 


(A. 9) 


Following  this  same  line  of  reasoning,  the  other  two  midplane  strain 

terms,  e°  and  7°  ,  can  be  written  as  [7] 
y  xy 

e  -  V,  +  I/2(u,  +  w,  +  V,  ) 

y  y  y  y  y 


(A. 10) 


7  -  U,  +  V,  +(u.  U,  +  V,  V,  +  w,  w,  ) 

xy  y  X  xy  xy  xy 


Now,  if  the  Kirchhoff  hypothesis  is  included,  where  strains  away  from 
the  midplane  are  due  to  the  following  curvatures 


K  — 

-w. 

X 

XX 

K  — 

-W, 

(A. 11) 

y 

yy 

K  — 

-w, 

xy 

xy 

the  full  expression  for  strain  in  the  plate  element  can  be  written 
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from  Equations  (A. 9),  (A. 10),  and  (A. 11)  as 


o 

e  -  €  -z  /c 

XX  X 


£ 

y 


7 


xy 


o 

€  -Z  K 

y  y 
7°  -2z  K 

xy  xy 


(A. 12) 


This  kinematic  formulation  allows  for  large  displacements  and 
moderate  rotations  (due  to  the  Kirchhoff  hypothesis) . 


Appendix  B:  Variational  Formulation 


STAGSC-1  is  an  energy  based  finite  element  program  dependent  on 
a  variational  formulation.  Therefore  a  presentation  of  this 
formulation  is  given  for  better  understanding. 

For  the  static  case  where  the  motion  of  the  system  can  be 
ignored,  the  total  potential  energy  of  a  system,  11,  can  be  given  as 

n  -  u  -  w  (B.i) 

where  U  is  the  internal  strain  energy  of  the  system  and  W  is  the 
external  work  due  to  applied  forces.  For  a  conservative  system,  the 

I 

change  in  the  total  potential  energy  is  independent  of  path  [7].  The 
equations  of  equilibrium  can  then  be  derived  from  the  first  variation 
of  the  total  potential  energy  [7].  In  equation  form,  the  principle 
of  virtual  work  is 

5n  -  5U  -  5W  -  0  (B.2) 


The  expression  for  the  internal  strain  energy  is  [6] 


-I 


(€)  [D]{€)  dVol 


(B.3) 


where  {«}  is  the  strain  vector  for  the  body  and  [D]  is  the  material 
matrix.  Taking  the  variation  of  Equation  (B.3)  results  in 


6U  -  1/2 


J  rt 


{€)  dVol  +  1/2 


V  t* 


]5{£)  dVol 


(B.4) 


or  alternatively  by  taking  the  transpose  of  the  last  expression  in 
Equation  (B.4)  and  adding  it  to  the  remaining  term  results  in 


5U 


lcr}dVol 


(B.5) 


Now  one  defines  a  differential  operator,  [L] ,  which  acts  on  the 
general  displacements,  (u) ,  to  get  [8] 
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le)  -  [L]{u}  (B.6) 

Next  one  can  define  the  shape  functions,  [N] ,  to  describe  the  general 
displacements  in  terms  of  nodal  displacements  (parameters),  (a),  of 
the  element  where 


{u)  -  [N]{a)  (B.7) 

If  one  substitutes  Equation  (B.7)  into  Equation  (B.6)  the  results 
give 

{€)  -  [L][N]{a)  (B.8) 

Or  alternatively,  defining  a  matrix  [B]  such  that  [8] 

[B]  -  (L1[N]  (B.9) 

Equation  (B.8)  can  now  be  written  as 


{€)  -  [B]{a}  (B.IO) 

Finally  one  can  take  the  variation  of  Equation  (B.IO),  transpose  it, 
and  substitute  it  into  Equation  (B.5)  resulting  in 

6U  -  «{a)^[  [B]^<7)  dVol  (B.ll) 

In  order  to  get  the  full  expression  for  the  variation  of  the 
potential  energy,  assume  that  the  work  term  has  been  defined  in  terms 
of  nodal  displacements,  (a),  and  forces,  {f),  such  that 

W  -  {a)^{f)  (B.12) 

Taking  the  variation  of  Equation  (B.12),  combining  it  with  Equation 
(B.ll),  and  substituting  these  equations  into  Equation  (B.2)  results 
in 


sn  -  Sla) 


1 H 


(o)  dVol  -  (f> 


}-0 


(B.13) 


Finally,  define  a  term,  {'t),  which  represents  the  sum  of  external  and 
internal  forces  (inside  braces.  Equation  (B.13))  [8] 
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(B.14) 


{*)  -  [B]^{a)  dVol  -  (f)  -  0 

■'v 

or  alternatively  after  integration 

{«)  -  [K]{a}  -  {f}  -  0  (B.15) 

where  [K]  is  the  stiffness  matrix. 
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Appendix  C:  Classical  Laminated  Plate  Theory 


The  shell  analyzed  in  this  thesis  is  made  of  composite 
materials.  STAGSC-1,  the  computer  program  used  in  the  analysis,  uses 
flat  plates  to  model  a  shell  surface.  Because  of  this,  a  review  of 
composite  plate  theory  is  presented. 

To  start  the  formulation,  a  few  of  the  classical  laminated  plate 
assumptions  need  to  be  made.  The  laminate  is  assumed  to  consist  of 
perfectly  bonded  plies  where  this  bond  is  assumed  infinitesimally 
chin  and  does  not  allow  shear  deformation  within  itself  (i.e.  the 
displacements  between  plies  is  continuous;  no  slip  between  plies) 
[14],  Also,  the  laminate  is  thin  so  that  the  in-plane  displacements, 
u  and  V,  as  well  as  the  normal  displacement,  w,  are  functions  of  only 
two  in-plane  space  variables  (x  and  y)  [7].  A  plane  stress 
assumption  is  assumed  to  where  7  ,  7  ,  c  ,  and  a  are  assumed  to 

xt  yz  z  z 

equal  zero.  Finally,  out  of  plane  strains  are  due  to  the  curvatures 
as  in  the  Kirchhoff  hypotheses  for  plates. 

Figure  C.l  shows  the  coordinate  system,  the  associated 
displacements,  and  the  force  and  moment  resultants  (N  ,  M  ,  etc. 

X  X 

respectively)  on  the  laminate  [14].  The  x-y  plane  of  Figure  C.l 
coincides  with  the  reference  axis.  From  these  assumptions,  a  general 
expression  can  be  written  for  the  strains  in  the  laminate  using  the 
Kirchhoff  hypothesis  as 
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where  €°,  e°,  and  7°  represent  in-plane  strains  and  the  terms  «  , 
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curvature  terms  are  defined  as 

(C.2) 


X  /  u 
Mxy 


Figure  C.l  Coordinate  System,  Displacements,  and 
Force  and  Moment  Resultants 


The  kinematic  relations  for  the  laminate  (Equation  C.l)  can  be 
substituted  into  the  constituitive  equation  for  the  laminate 
resulting  in  [1^] 
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^  xy^ 

* 

(C.3) 


where  the  subscript  k  refers  to  the  ply  layer  and  the  Q  quantities 
are  defined  as  follows  [1^]: 


-  Q  cos^fl  +  2(Q  +  2Q  )sin^^  +  Q  sin*^ 

11  11  12  ^66'  ^22 

Q  -  (Q  +  Q  -  4Q  )sin^^  cos^tf  +  Q  (sin*S  +  cos*^) 
12  11  22  66  ^12 


(C.4) 
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(C.4) 

(cont) 


Q  -  Q  sin*^  +  2(Q  +  2Q  )sln^5  cos^B  +  Q  cos*$ 

22  11  12  66  22 

Q  “  (Q  *  Q  -  2Q  )sin  6  cos^ B  +  (Q  -  Q  +  Q  )sin^S  cos  B 

^16  ^12  ^66"^  '^12  ^22  ^66 

Q  -  (Q  -  Q  -  2Q  )sin^d  cos  5  +  (Q  -  Q  +  Q  )sin  B  cos^ B 

26  II  12  66  12  22  66 

Q  -  (Q  -  Q  -  2Q  -  2Q  )sin^fl  cos^$  +  Q  (sin*^  +  cos*^) 

66  11  22  16  6  6  66 


The  angle  B  in  Equation  (C.4)  is  the  angle  between  the  principle  axis 
(x,y,z)  and  the  material  axis  (1,2,3)  as  in  Figure  C.2 

The  force  and  moment  resultants  (forces  and  moments  per  unit 
length)  acting  on  a  composite  plate  can  be  found  by  integrating  the 
stresses  in  each  ply  through  the  plate  thickness,  for  example  [14] 

t*2 

N  -  la  dz 

X  I  X 

-t7  2 

(C.6) 

t|.2  ' 

M  -  I  a  z  dz 

X  J  X 

-t72 

These  force  and  moment  resultants  are  shown  in  Figure  C.l.  All  of 
the  force  and  moment  resultants  are  put  in  vector  form  for  an 
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Figure  C.2  Principal  (x,y,z)  and  Material  (1,2,3) 
Axis  Systems 


N- layered  laminate  and  the  results  are  [lA] 
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where  z  and  z  are  defined  in  Figure  C.3. 
k  k-l  ^ 


(C.7) 


(C.8) 


Substituting  Equation  (C.3)  into  Equations  (C.7)  and  (C.8),  results 
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Figure  C.3  Geometry  of  an  N-Layered  Laminate  [14] 

The  stiffness  matrix  [Q]  has  been  removed  from  the  Integrals  in 
Equations  (C.9)  because  it  is  constant  within  a  lamina  but  must  be 
included  within  the  summation  of  force  and  moment  resultants  for  each 
ply  [14],  Since  the  strain  and  curvature  terms  in  Equations  (C.9) 
are  not  functions  of  z  but  are  midplane  values  they  can  be  removed 
from  the  integration  and  summation  signs  [14].  After  integrating  the 
remaining  terms,  the  results  are 
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Appendix  D:  Surface  Normal  Calculations 


In  order  to  compare  the  global  displacements  calculated  by  STAGS 
with  the  normal  displacements  from  experimentation,  the  STAGS 
displacements  had  to  be  transformed  into  surface  normal 
displacements.  To  do  this  an  approximate  unit  normal  was  calculated 
at  each  node  point  on  the  finite  element  model  and  multiplied  (scalar 
product)  by  the  global  displacements  (in  vector  form)  resulting  in  a 
normal  displacement. 

An  example  of  a  unit  normal  calculation  for  a  hypothetical  node 
(node  5)  in  Figure  D.l  is  shown. 


Figure  D.l  Hypothetical  Finite  Element  Surface 

First,  the  element  dividing  lines  radiating  from  node  5  to  the 
adjacent  nodes  are  converted  into  vectors.  For  example,  to  get  a 
vector  from  node  5  to  node  1  (b  )  the  global  coordinates  of  node  1 


(g^)  are  subtracted  as  a  vector  from  the  global  coordinates  of  node  5 
(g^)  as  in 


Next,  the  vectors  radiating  out  from  node  5  are  multiplied 
(cross  product)  for  each  element  (numbers  circled  in  Figure  D.l)  to 
find  a  normal  for  the  element.  For  example,  for  element  4 

a-bxb  (D.2) 

41  54  51  ^  ' 

where  the  resulting  vector  a^^  is  an  outward  normal. 

The  four  vectors  (a  ,  a  ,a  ,  and  a  )  are  then  normalized, 
added  together,  and  divided  by  four  resulting  in  a  normal  vector 
radiating  from  the  node.  Finally,  this  vector  is  then  normalized  and 
multiplied  (scalar  product)  by  the  global  displacements  (in  vector 
form)  resulting  in  an  approximate  normal  displacement. 

A  Fortran  program  was  written  to  calculate  the  surface  unit 
normals  for  the  finite  element  model  and  stored  in  a  file.  As  the 
displacement  solutions  were  available,  they,  as  wel],  as  the  unit 
normals,  were  used  as  input  files  to  another  Fortran  program  that 
calculated  normal  displacements  for  the  entire  shell. 
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Appendix  E:  Convergence  Model;  1  to  1  Aspect  Ratio 


The  following  input  data  deck  was  used  for  the  1  to  1  aspect 
ratio  convergence  model.  The  reader  is  directed  to  the  STAGS  users 
manual  [16]  for  specifics  on  the  data  deck.  The  letters  immediately 
following  the  dollar  sign  in  each  row  can  be  cross  referenced  to 
specific  cards  in  the  STAGS  users  manual. 


(CONST  THICK,  COMPOSITE) 

$B1  ANALYSIS  TYPE 
SHELL  UNITS 
DATA  TABLES  SUMMARY 
0  $C1  LOAD  INCREMENT 

STRATEGY  PARAMETER  RCD 


NON -LIN  ANAL  CONVERGl-1 
3  1110  0  1 
1 

10  1 

0.1  0.1  1.0  0.0  0.0  0.0  0 
0  1800  2  -20  1 
37  15 
1 

4.89E06  0.035  1.9E05  0.0522  1.0  4.23E05  1.0  $12  MATERIAL  PROP 


$B2 
$B3 
0  0 
$D1 
$F1 

$11  KEVLAR/POLY  MATERIAL  RECORD 


1 

1  40 

$K1 

SHELL  WALL  RECORD 

1 

0.0018000 

0.0  1 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0018000 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0018000 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0018000 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0018000 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0018000 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0018000 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0018000 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 
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1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0018000 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0018000 

90.0  1 

$K2 

LAYERED 

WALL 

RECORD 

11 

0. 

1 

410 
0. 
27. 
152 
172 
2 
2 
2 
2 
1 
1 
0 


0 

180. 


0.  180.  108.  26. 


-1 

2.87  5.73  7.5  9.26  10 
81  33.71  39.6  49.8  60 
.19  155.52  158.84  161 
.5  174.27  177.13  180. 


0  2 


$M1 

$M2A 

$M5 

$N1 

,98  12.7  14.42  16.14  18. 
,  75.  90.  105.  120.  130. 
.35  163.86  165.58  167.3 
$N4 

$P1 

$P1 

$P1 

$P1 

$Q1 

$Q2 

$R1 


65  21.16  24.5, 
2  140.4  146.3, 
169.02  170.7, 
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Appendix  F:  UFRESS  Subroutine  for  Convergence  Test 

The  following  subroutine  was  used  for  applying  the  pressure 
loading  in  the  convergence  test.  The  subroutine  is  written  using 
guidelines  in  the  STAGS  users  manual  [16].  The  subroutine  is  then 
compiled  and  linked  to  the  main  program  before  execution. 


subroutine  upress(t ,pa,pb, iunit, ielt ,x,y , z , live, press) 
c 

c  Pressure  distribution  subroutine  for 
c  convergence  test.  Pressures  are  fed  back  to  the  main 
c  program  depending  on  element  x  and  y  coordinates . 
c  loads.dat  is  the  pressure  distribution  data  file 
c  that  has  the  twenty  loads,  one  on  each  line 
c 

integer  i 
real  c(20) 

open(unit-l , file-' loads . dat' , status-' old' ) 
live-1 

do  20  i-1.20 
read<l,30)  c(i) 

20  continue 

30  forraat(fl5.7) 

rewind (unit-1) 
c 

if  (x.le.16.14)  then 
if  (y.le.90.)  then 
press-c(l)*pa 
elseif  (y.le.180.)  then 
press-c(ll)*pa 
endif 

elseif  (x.le.39.6)  then 
if  (y.le.45.)  then 
press-c(3)*pa 
elseif  (y.le.90.)  then 
press-c(2)*pa 
elseif  (y.le.l35.)  then 
press-c(12)*pa 
elseif  (y.le.180.)  then 
press-c(13)*pa 
endif 

elseif  (x.le.90.)  then 
if  (y.le.45.)  then 
press-c(5)*pa 
elseif  (y.le.90.)  then 
press-c (4)*pa 
elseif  (y.le.l35.)  then 
press-c ( 14)*pa 
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elself  (y.le.180.)  then 
press-c(15)*pa 
endif 

elseif  (x.le.140.4)  then 
if  (y.le.45.)  then 
press-c(7)*pa 
elseif  (y.le.90.)  then 
press-c(6)*pa 
elseif  (y.le.l35.)  then 
press-c(16)*pa 
elseif  (y.le.l80.)  then 
press-c(17)*pa 
endif 

elseif  (x. le . 163 . 86)  then 
if  (y.le.45.)  then 
press-c(9)*pa 
elseif  (y.le.90.)  then 
press-c(8)*pa 
elseif  (y.le.l35.)  then 
press-c(18)*pa 
elseif  (y.le.180.)  then 
press-c(19)*pa 
endif 

elseif  (x.le.l80.)  then 
if  (y.le.90.)  then 
press-c(10)*pa 
elseif  (y.le.180.)  then 
press-c(20)*pa 
endif 
endif 
return 
endg 
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Appendix  G;  Finite  Element  Model  (Model  6) 


The  following  input  data  deck  was  used  for  the  Model  6  (see 
section  3.3)  finite  element  rtin.  The  reader  is  directed  to  the  STAGS 
users  manual  [16]  for  specifics  on  the  data  deck.  The  letters 
immediately  following  the  dollar  sign  in  each  row  can  be  cross 
referenced  to  specific  cards  in  the  STAGS  users  manual. 


NON-LIN  ANAL  UPRESS  LOAD  (VAR  THICK,  COMPOSITE,  NEW  B.C.) 
3111001  $B1  ANALYSIS  TYPE 

0  1  $B2  SHELL  UNITS 

41270  $B3  DATA  TABLES  SUMMARY 

0.1  0.1  1.0  0.0  0.0  0.0  0  0  0  0  $C1  LOAD  INCREMENT 

0  1800  2  -20  1  $D1  STRATEGY  PARAMETER  RCD 

466  0  52  54  360  0  0  $H1  ELEMENT  UNIT  SUMMARY 

1  $11  KEVLAR/POLY  MAT  RECORD 
4.89E06  0.035  1.9E05  0.0522  1.0  4.23E05  1.0  $12  MAT  PROP 

2  $11  GLASS  MAT  RECORD  (ISOTROPIC) 

3.24E06  0.33  0.0  0.050608  1.0  3.24E06  1.0  $12  MAT  PROP 

3  $11  ALUMINUM  MAT  RECORD  (ISOTROPIC) 
10.7E06  0.33  0.0  0.1  1.0  10.7E06  1.0  $12  MAT  PROP 

4  $12  RUBBER  BEAM  MAT  RECORD 

2.071E04  0.5  0.0  0.01  1.0  2.071E04  1.0  $12  MAT  PROP 
114000. 00. 00  $J1  CROSS  SECTION  RECORD 

1.146E-01  1.045E-03  1.045E-03  0.0  $J2A  CROSS  SECTION  RECORD 


1 

1  32 

$K1 

SHELL  WALL  RECORD 

1 

0.0042500 

0.0  1 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

0.0 

$K2 

LAYERED 

WALL 

RECORD 

1 

0.0042500 

90.0 

$K2 

LAYERED 

WALL 

RECORD 
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1  0.0042500 

0.0 

$K2  LAYERED  WALL  RECORD 

1  0.0042500  90.0 

$K2  LAYERED  WALL  RECORD 

1  0.0042500 

0.0 

$K2  LAYERED  WALL  RECORD 

1  0.0042500  90.0 

$K2  LAYERED  WALL  RECORD 

1  0.0042500 

0.0 

$K2  LAYERED  WALL  RECORD 

1  0.0042500  90.0 

$K2  LAYERED  WALL  RECORD 

1  0.0042500 

0.0 

$K2  LAYERED  WALL  RECORD 

1  0.0042500  90.0 

$K2  LAYERED  WALL  RECORD 
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3 

0 

0 

0 

1.7067 

2.6154 

2.4038 

111 

111  $S1  USER 

PT 
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4 

0 

0 

0 

1.7067 

.0000 

3.4593 

111 

111 

$S1  USER 

PT 

RCD 

5 

0 

0 

0 

1.7067 

2.6154 

-2.4038 

111 

111 
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PT 

RCD 

6 

0 

0 

0 

1.7067 

.0000 

-3.4593 

111 

111  $S1  USER 

PT 

RCD 
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0 

0 

0 

5.8027 

9.2129 
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111 

111 
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PT 

RCD 

8 

0 

0 

0 

5.8027 

8.0588 

2.7373 

111 

111  $S1  USER 

PT 
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9 

0 

0 

0 

5.8027 

5.8590 

4.8204 

111 

111 
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PT 

RCD 

10 

0 

0 

0 

5.8027 

2.4660 

6.7054 

111 

111 

$S1  USER 

PT 

RCD 

11 

0 

0 

0 

5.8027 

-  .1554 

6.9066 

111 

111 

$S1  USER 

PT 

RCD 

12 

0 

0 

0 

5.8027 

8.0588 

-2.7373 

111 

111 

$S1  USER 

PT 

RCD 

13 

0 

0 

0 

5.8027 

5.8590 

-4.8204 

111 

111 

$S1  USER 

PT 

RCD 

14 

0 

0 

0 

5.8027 

2.4660 

-6.7054 

111 

111 
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PT 

RCD 

15 

0 

0 

0 

5.8027 

- .1554 

-6.9066 

111 

111 
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PT 

RCD 
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0 

0 

0 

11.6267 

13.1878 
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111 

111 
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PT 
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17 

0 

0 

0 

11.6267 

12.5106 

2.7447 

111 

111 
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PT 
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18 

0 

0 

0 

11.6267 

10  9231 

5.1033 

111 

111 
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PT 

RCD 

19 

0 

0 

0 

11.6267 

8.8064 

7.0056 

111 

111 

$S1  USER 

PT 

RCD 

20 

0 

0 

0 

11.6267 

6.3305 

8.4114 

111 

111 

$S1  USER 

PT 

RCD 

21 

0 

0 

0 

11.6267 

2.2535 

9.4921 

111 

111 

$S1  USER 

PT 

RCD 

22 

0 

0 

0 

11.6267 

- .3141 

9.8173 

111 

111 

$S1  USER 

PT 

RCD 

23 

0 

0 

0 

11.6267 

12.5106 

-2.7447 

111 

111 

$S1  USER 

PT 

RCD 

24 

0 

0 

0 

11.6267 

10.9231 

-5.1033 

111 

111 

$S1  USER 

PT 

RCD 

25 

0 

0 

0 

11.6267 

8.8064 

-7.0056 

111 

111 

$S1  USER 

PT 

RCD 

26 

0 

0 

0 

11.6267 

6.3305 

-8.4114 

111 

111 

$S1  USER 

PT 

RCD 

27 

0 

0 

0 

11.6267 

2.2535 

-9.4921 

111 

111 

$S1  USER 

PT 

RCD 

28 

0 

0 

0 

11.6267 

- .3141 

-9.8173 
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111 
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PT 

RCD 

29 

0 

0 

0 

17.4507 

16.4566 

.0000 

111 

111 

$S1  USER 

PT 

RCD 

30 

0 

0 

0 

17.4507 

16.1636 

2.2341 

111 

111 

$S1  USER 

PT 

RCD 

31 

0 

0 

0 

17.4507 

15.3482 

4.3357 

111 

Hi 

$S1  USER 

PT 

RCD 

32 

0 

0 

0 

17.4507 

13.2931 

7.0699 

111 

111 

$S1  USER 

PT 

RCD 

33 

0 

0 

0 

17.4507 

10.5753 

9.1282 

111 

111 

$S1  USER 

PT 

RCD 

34 

0 

0 

0 

17.4507 

7.4389 

10.5816 

111 

111 

$S1  USER 

PT 

RCD 

35 

0 

0 

0 

17.4507 

2.0410 

11.6746 

111 

111 

$S1  USER 

PT 

RCD 

36 

0 

0 

0 

17.4507 

- .4632 

11.9202 

111 

111 

$S1  USER 

PT 

RCD 

37 

0 

0 

0 

17.4507 

16.1636 

-2.2341 

111 

111 

$S1  USER 

PT 

RCD 

38 

0 

0 

0 

17.4507 

15.3482 

-4.3357 

111 

111 

$S1  USER 

PT 

RCD 

39 

0 

0 

0 

17.4507 

13.2931 

-7.0699 

111 

111 
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PT 

RCD 

40 

0 

0 

0 

17.4507 

10.5753 

-9.1282 

111 
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PT 

RCD 

41 

0 

0 
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17.4507 

7.4389 

-10.5816 
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$S1  USER 
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0 
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17.4507 

2.0410 

-11.6746 
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PT 
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0 

0 
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- .4632 

-11.9202 
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0 
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23.2533 

19.3166 
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45 
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0 
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23.2533 

18.9387 

2.9493 

111 

111 
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PT 

RCD 

46 

0 

0 

0 

23.2533 

17.7758 

5.6968 

111 

111 
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PT 

RCD 

47 

0 

0 

0 

23.2533 

15.1888 

8.6430 

111 

111 
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PT 

RCD 

48 

0 

0 

0 

23.2533 

11.9307 

10.7546 
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PT 
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49 
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0 
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23.2533 

8.2676 

12.2041 
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4.1848 

13.0923 
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0 

0 

0 

23.2533 

1.9000 

13.4892 
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PT 
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0 

0 

0 

23.2533 

-.5935 

13.4892 

111 
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$S1 
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PT 

ROD 

53 

0 

0 

0 

23.2533 

18.9387 

-2.9493 

111 
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$S1 
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PT 

RCD 

54 

0 

0 

0 

23.2533 

17.7758 

-5.6968 

111 
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$S1 
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PT 

RCD 

55 

0 

0 
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23.2533 

15.1888 

-8.6430 
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$S1 
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PT 
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23.2533 

11.9307 

-10.7546 
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4.1848 
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14,2577 
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9.3745 
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4.4698 

15.0336 
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23.4369 
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21.6871 
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14.2577 
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9.3745 

-14.3713 
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4.4698 

-15.0336 
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1.6900 

-15.7464 
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46.5493 
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46.5493 

26.6173 

4.9282 
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46.5493 

24.4760 

9.4544 
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111 

$S1 
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81 
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0 
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46.5493 

20.7198 

13.2326 
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46.5493 

15.7233 
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46.5493 

10.2022 
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46.5493 

4.7007 
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46.5493 

1.5300 
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46.5493 
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46.5493 
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-9.4544 
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46.5493 

20.7198 
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46.5493 
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46.5493 

4.7007 
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46.5493 
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58.1760 
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28.7711 

5.4832 

111 

111 

$S1 
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97 

0 

0 

0 
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26.3311 

10.4938 
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98 

0 

0 

0 

58.1859 

21.9052 

14.2165 
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RCD 

99 

0 

0 

0 

58.1759 

16.4688 

16.2003 

111 
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$S1 
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PT 
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100 
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0 

0 

58.1759 

10.6657 

16.8997 

111 

111 

$S1 
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58.1759 

4.8376 

17.5387 
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1.4200 
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18.1876 

111 

111 

$S1 

USER 

PT 

RCD 
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58.1759 

28.7711 
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58.1759 

26.3311 

-10.4938 
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$S1 
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58.1759 
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-14.2165 
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58.1759 

16.4688 
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10.6657 
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58.1759 

4.8376 
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110 

0 
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58.1759 

1.4200 

-18.1876 
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0 
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58.1759 

-1.0806 
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69.8240 
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30.2019 

5.8520 
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0 

0 

0 

69.8239 
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11.0656 
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0 
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69.8239 

16 . 9340 

16.7585 
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0 
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17.4168 
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$S1 

USER 

PT 

RCD 
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0 

0 

0 

69.8239 

4.9165 

18.0761 

111 
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119 

0 

0 

0 

69.8239 

1.3500 

18 . 7409 

111 
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PT 
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0 

0 

0 

69.8239 

-1.1476 

18.7409 
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$S1 
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PT 
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0 

0 

0 

69.8239 

30.2019 

-5.8520 

111 
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$S1 
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PT 
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0 

0 

0 

69.8239 

27.3508 

-11.0656 
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$S1 
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0 

69.8239 

22.6940 

-14.8710 
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69.8239 

16.9340 

-16.7585 
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10.9299 

-17.4168 

111 

111 

$S1 

USER 

PT 

RCD 

126 

0 

0 

0 

69.8239 

4.9165 

-18.0761 
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69.8239 
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-18.7409 

111 

111 

$S1 

USER 

PT 

RCD 

128 

0 

0 

0 
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75.6273 
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111 

111 

$S1 
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0 

75.6373 

30.5637 

5.9453 

111 
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$S1 
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PT 
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0 

0 

0 

75.6373 

27.5961 

11.2032 

111 

111 

$S1 

USER 

PT 

RCD 

132 

0 

0 

0 

75.6373 

22.8776 

15.0234 

111 

111 

$S1 

USER 

PT 

RCD 

133 

0 

0 

0 

75.6373 

17.0874 

16.9425 

111 

111 

$S1 

USER 

PT 

RCD 

134 

0 

0 

0 

75.6373 

11.0153 

17.5839 

111 

111 

$S1 
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PT 

RCD 

135 

0 

0 

0 

75.6373 

4.9385 

18.2258 

111 

111 

$S1 
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PT 

RCD 
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0 

0 

75.6373 

1 . 3400 

18.8704 
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75.6373 

-1.1636 

18.8704 
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17.0874 
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11.0153 

-17.5839 
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PT 
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75.6373 

4.9385 

-18.2258 
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PT 
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0 

75.6373 

1.3400 

-18.8704 
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$S1 
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PT 
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0 
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0 

75.6373 

-1.1636 

-18.8704 
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PT 
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0 

0 

81.4506 

31.9999 
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0 

81.4506 

30.9256 

6.0385 
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PT 
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0 

0 

81.4506 

27.8414 

11.3408 
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PT 
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81.4506 

23.0612 

15.1759 
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$S1 

USER 

PT 

RCD 
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81.4506 

17.2408 

17.1266 
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PT 
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81.4506 

11.1006 

17.7510 
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PT 
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81.4506 

4.9605 

18.3755 
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111 

111 

$S1 

USER 

PT 

ROD 

159 

0 

0 

0 

81.4506 

11.1006  -17.7510 

111 

111 

$S1 

USER 

PT 

RCD 

160 

0 

0 

0 

81.4506 

4.9605  -18.3755 

111 

111 

$S1 

USER 

PT 

RCD 

161 

0 

0 

0 

81.4506 

1.3200  -18.9999 

111 

111 

$S1 

USER 

PT 

RCD 

162 

0 

0 

0 

81.4506 

-1.1797  -18.9999 

111 

111 

$S1 

USER 

PT 

RCD 

163 

0 

0 

0 

90.8613 

32.0000 

.0000 

111 

111 

$S1 

USER 

PT 

RCD 

164 

0 

0 

0 

90.8613 

30.8852 

6.0320 

111 

111 

$S1 

USER 

PT 

RCD 

165 

0 

0 

0 

90.8613 

27.7692 

11.3113 

111 

111 

$S1 

USER 

PT 

RCD 

166 

0 

0 

0 

90.8613 

23.0409 

15.1589 

111 

111 

$S1 

USER 

PT 

RCD 

167 

0 

0 

0 

90.8613 

17.2351 

17.1560 

111 

111 

$S1 

USER 

PT 

RCD 

168 

0 

0 

0 

90.8613 

11.1349 

17.7957 

111 

111 

$S1 

USER 

PT 

RCD 

169 

0 

0 

0 

90.8613 

4.9926 

18.4399 

111 

111 

$S1 

USER 

PT 

RCD 

170 

0 

0 

0 

90.8613 

1.3100 

18.9331 

111 

111 

$S1 

USER 

PT 

RCD 

171 

0 

0 

0 

90.8613 

-1.1880 

19.0881 

111 

111 

$S1 

USER 

PT 

RCD 

172 

0 

0 

0 

90.8613 

30.8852 

-6.0320 

111 

111 

$S1 

USER 

PT 

RCD 

173 

0 

0 

0 

90.8613 

27.7692  -11.3113 

111 

111 

$S1 

USER 

PT 

RCD 

174 

0 

0 

0 

90.8613 

23.0409  -15.1589 

111 

111 

$S1 

USER 

PT 

RCD 

175 

0 

0 

0 

90.8613 

17.2351  -17.1560 

111 

111 

$S1 

USER 

PT 

RCD 

176 

0 

0 

0 

90.8613 

11.1349  -17.7957 

111 

111 

$S1 

USER 

PT 

RCD 

177 

0 

0 

0 

90.8613 

4.9926  -18.4399 

111 

111 

$S1 

USER 

PT 

RCD 

178 

0 

0 

0 

90.8613 

1.3100  -18.9331 

111 

111 

$S1 

USER 

PT 

RCD 

179 

0 

0 

0 

90.8613 

-1.1880  -19.0881 

111 

111 

$S1 

USER 

PT 

RCD 

180 

0 

0 

0 

98.5C00 

32.0000 

.0000 

111 

111 

$S1 

USER 

PT 

RCD 

181 

0 

0 

0 

98.5000 

30.8852 

6.0320 

111 

111 

$S1 

USER 

PT 

RCD 

182 

0 

0 

0 

98.5000 

27.7692 

11.3113 

111 

111 

$S1 

USER 

PT 

RCD 

183 

0 

0 

0 

98.5000 

23.0409 

15.1589 

111 

111 

$S1 

USER 

PT 

RCD 

184 

0 

0 

0 

98.5000 

17.2351 

17.1560 

111 

111 

$S1 

USER 

PT 

RCD 

185 

0 

0 

0 

98.5000 

11.1349 

17.7957 

111 

111 

$S1 

USER 

PT 

RCD 

186 

0 

0 

0 

98.5000 

4.9926 

18.4399 

111 

111 

$S1 

USER 

PT 

RCD 

187 

0 

0 

0 

98.5000 

1.1300 

18.9331 

111 

111 

$S1 

USER 

PT 

RCD 

188 

0 

0 

0 

98.5000 

-1.1880 

19.0881 

111 

111 

$S1 

USER 

PT 

RCD 

189 

0 

0 

0 

98.5000 

30.8852 

-6.0320 

111 

111 

$S1 

USER 

PT 

RCD 

190 

0 

0 

0 

98.5000 

27.7692  -11.3113 

111 

111 

$S1 

USER 

PT 

RCD 

191 

0 

0 

0 

98 . 5000 

23.0409  -15.1589 

111 

111 

$S1 

USER 

PT 

RCD 

192 

0 

0 

0 

98.5000 

17.2351  -17.1560 

111 

111 

$S1 

USER 

PT 

RCD 

193 

0 

0 

0 

98.5000 

11.1349  -17.7957 

111 

111 

$S1 

USER 

PT 

RCD 

194 

0 

0 

0 

98.5000 

4.9926  -18.4399 

111 

111 

$S1 

USER 

PT 

RCD 

195 

0 

0 

0 

98.5000 

1.3100  -18.9331 

111 

111 

$S1 

USER 

PT 

RCD 

196 

0 

0 

0 

98.5000 

-1.1880  -19.0881 

111 

111 

$S1 

USER 

PT 

RCD 

197 

0 

0 

0 

103.8900 

31.9687 

.0000 

111 

111 

$S1 

USER 

PT 

RCD 

198 

0 

0 

0 

103.8900 

30.8852 

6.0320 

111 

111 

$S1 

USER 

PT 

RCD 

199 

0 

0 

0 

103.8900 

27.7692 

11.3113 

111 

111 

$S1 

USER 

PT 

RCD 

200 

0 

0 

0 

103.8900 

23.0409 

15.1589 

111 

111 

$S1 

USER 

PT 

RCD 

201 

0 

0 

0 

103.8900 

17.2351 

17.1560 

111 

111 

$S1 

USER 

PT 

RCD 

202 

0 

0 

0 

103.8900 

11.1349 

17.7957 

111 

111 

$S1 

USER 

PT 

RCD 

203 

0 

0 

0 

103.8900 

4.9926 

18.4399 

111 

111 

$S1 

USER 

PT 

RCD 

204 

0 

0 

0 

103.8900 

1.3100 

18.9331 

111 

111 

$S1 

USER 

PT 

RCD 

205 

0 

0 

0 

103.8900 

-1.1880 

19.0881 

111 

111 

$S1 

USER 

PT 

RCD 

206 

0 

0 

0 

103.8900 

30.8852 

-6.0320 

111 

111 

$S1 

USER 

PT 

RCD 

207 

0 

0 

0 

103.8900 

27.7692  -11.3113 

111 

111 

$S1 

USER 

PT 

RCD 

208 

0 

0 

0 

103.8900 

23.0409  -15.1589 

111 

111 

$S1 

USER 

PT 

RCD 

209 

0 

0 

0 

103.8900 

17.2351  -17.1560 

111 

111 

$S1 

USER 

PT 

RCD 

210 

0 

0 

0 

103.8900 

11.1349  -17.7957 

111 

111 

$S1 

USER 

PT 

RCD 

211 

0 

0 

0 

103.8900 

4.9926  -18.4399 

111 

111 

$S1 

USER 

PT 

RCD 

141 


212 

0 

0 

0 

103.8900 

1.3100 

-18.9331 

111 

111 

$S1 

USER 

PT 

ROD 

213 

0 

0 

0 

103.8900 

-1.1880 

-19.0881 

111 

111 

$S1 

USER 

PT 

ROD 

214 

0 

0 

0 

113.3800 

31.4000 

.0000 

111 

111 

$S1 

USER 

PT 

ROD 

215 

0 

0 

0 

113.3800 

30.2179 

5.9301 

111 

111 

$S1 

USER 

PT 

ROD 

216 

0 

0 

0 

113.3800 

27.2238 

11.0903 

111 

111 

$S1 

USER 

PT 

ROD 

217 

0 

0 

0 

113.3800 

22.7083 

14.8830 

111 

111 

$S1 

USER 

PT 

ROD 

218 

0 

0 

0 

113.3800 

17.0730 

16.9586 

111 

111 

$S1 

USER 

PT 

ROD 

219 

0 

0 

0 

113.3800 

11.0068 

17.6280 

111 

111 

$S1 

USER 

PT 

ROD 

220 

0 

0 

0 

113.3800 

4.9016 

18.2567 

111 

111 

$S1 

USER 

PT 

RCD 

221 

0 

0 

0 

113.3800 

1.3300 

18.7300 

111 

111 

$S1 

USER 

PT 

ROD 

222 

0 

0 

0 

113.3800 

-1.1686 

18.8819 

111 

111 

$S1 

USER 

PT 

RCD 

223 

0 

0 

0 

113.3800 

30.2179 

-5.9301 

111 

111 

$S1 

USER 

PT 

RCD 

224 

0 

0 

0 

113.3800 

27.2238 

-11.0904 

111 

111 

$S1 

USER 

PT 

RCD 

225 

0 

0 

0 

113.3800 

22.7083 

-14.8831 

111 

111 

$S1 

USER 

PT 

RCD 

226 

0 

0 

0 

113.3800 

17.0730 

-16.9586 

111 

111 

$S1 

USER 

PT 

RCD 

227 

0 

0 

0 

113.3800 

11.0068 

-17.6280 

111 

111 

$S1 

USER 

PT 

RCD 

228 

0 

0 

0 

113.3800 

4.9016 

-18.2567 

111 

111 

$S1 

USER 

PT 

RCD 

229 

0 

0 

0 

113.3800 

1.3300 

-18.7300 

111 

111 

$S1 

USER 

PT 

RCD 

230 

0 

0 

0 

113.3800 

-1.1686 

-18.8819 

111 

111 

$S1 

USER 

PT 

RCD 

231 

0 

0 

0 

122.8700 

30.5737 

.0000 

111 

111 

$S1 

USER 

PT 

RCD 

232 

0 

0 

0 

122.8700 

29.5506 

5.8281 

111 

111 

$S1 

USER 

PT 

RCD 

233 

0 

0 

0 

122.8700 

26.6784 

10.8694 

111 

111 

$S1 

USER 

PT 

RCD 

234 

0 

0 

0 

122.8700 

22.3758 

14.6072 

111 

111 

$S1 

USER 

PT 

RCD 

235 

C 

0 

0 

122.8700 

16.9109 

16.7612 

111 

111 

$S1 

USER 

PT 

RCD 

236 

0 

0 

0 

122.8700 

10.8788 

17.4604 

111 

111 

$S1 

USER 

PT 

RCD 

237 

0 

0 

0 

122.8700 

4.8107 

18.0735 

111 

111 

$S1 

USER 

PT 

RCD 

238 

0 

0 

0 

122.8700 

1.-500 

18.5269 

111 

111 

$S1 

USER 

PT 

RCD 

239 

0 

0 

0 

122.8700 

-1.1492 

18.6758 

111 

111 

$S1 

USER 

PT 

RCD 

240 

0 

0 

0 

122.8700 

29.5506 

-5.8281 

111 

111 

$S1 

USER 

PT 

RCD 

241 

0 

0 

0 

122.8700 

26.6784 

-10.8694 

111 

111 

$S1 

USER 

PT 

RCD 

242 

0 

0 

0 

122.8700 

22.3758 

-14.6072 

111 

111 

$S1 

USER 

PT 

RCD 

243 

0 

0 

0 

122.8700 

16.9109 

-16.7612 

111 

111 

$S1 

USER 

PT 

RCD 

244 

0 

0 

0 

122.8700 

10.8788 

-17.4604 

111 

111 

$S1 

USER 

PT 

RCD 

245 

0 

0 

0 

122.8700 

4.8107 

-18.0735 

111 

111 

$S1 

USER 

PT 

RCD 

246 

0 

0 

0 

122.8700 

1.3500 

-18.5269 

111 

111 

$S1 

USER 

PT 

RCD 

247 

0 

0 

0 

122.8700 

-1.1492 

-18.6758 

111 

111 

$S1 

USER 

PT 

RCD 

248 

0 

0 

0 

136.1700 

28.2630 

.0000 

111 

111 

$S1 

USER 

PT 

RCD 

249 

0 

0 

0 

136.1700 

27.3259 

5.4883 

111 

111 

$S1 

USER 

PT 

RCD 

250 

0 

0 

0 

136.1700 

24.8112 

10.1130 

111 

111 

$S1 

USER 

PT 

RCD 

251 

0 

0 

0 

136.1700 

21.1503 

13.5906 

111 

111 

$S1 

USER 

PT 

RCD 

252 

0 

0 

0 

136.1700 

16.2183 

15.9177 

111 

111 

$S1 

USER 

PT 

RCD 

253 

0 

0 

0 

136.1700 

10.3390 

16.7532 

111 

111 

$S1 

USER 

PT 

RCD 

254 

0 

0 

0 

136.1700 

4.4403 

17.3274 

111 

111 

$S1 

USER 

PT 

RCD 

255 

0 

0 

0 

136.1700 

1.4500 

17.7217 

111 

111 

$S1 

USER 

PT 

RCD 

256 

0 

0 

0 

136.1740 

-1.0464 

17.8586 

111 

111 

$S1 

USER 

PT 

RCD 

257 

0 

0 

0 

136.1700 

27.3259 

-5.4883 

111 

111 

$S1 

USER 

PT 

RCD 

258 

0 

0 

0 

136.1700 

24.8112 

-10.1130 

111 

111 

$S1 

USER 

PT 

RCD 

259 

0 

0 

0 

136.1700 

21.1503 

-13.5906 

111 

111 

$S1 

USER 

PT 

RCD 

260 

0 

0 

0 

136.1700 

16.2183 

-15.9177 

111 

111 

$S1 

USER 

PT 

RCD 

261 

0 

0 

0 

136.1700 

10.3390 

-16.7532 

111 

111 

$S1 

USER 

PT 

RCD 

262 

0 

0 

0 

136.1700 

4.4403 

-17.3274 

111 

111 

$S1 

USER 

PT 

RCD 

263 

0 

0 

0 

136.1700 

1.4500 

-17.7217 

111 

111 

$S1 

USER 

PT 

RCD 

264 

0 

0 

0 

136.1740 

-1.0464 

-17.8586 

111 

111 

$S1 

USER 

PT 

RCD 

265 

0 

0 

0 

149.4800 

24.9563 

.0000 

111 

111 

$S1 

USER 

PT 

RCD 

142 


266 

0 

0 

0 

149.4800 

24.1245 

4.9995 

111 

111 

$S1 

USER 

PT 

RCD 

267 

0 

0 

0 

149.4800 

22.0625 

8.9998 

111 

111 

$S1 

USER 

PT 

RCD 

268 

0 

0 

0 

149.4800 

19.2328 

11.9998 

111 

111 

$S1 

USER 

PT 

RCD 

269 

0 

0 

0 

149.4800 

15.0000 

14.4343 

111 

111 

$S1 

USER 

PT 

RCD 

270 

0 

0 

0 

149.4800 

9.3972 

15.5198 

111 

111 

$S1 

USER 

PT 

RCD 

271 

0 

0 

0 

149.4800 

3.8030 

16.0437 

111 

111 

$S1 

USER 

PT 

RCD 

272 

0 

0 

0 

149.4800 

1.5600 

16.3084 

111 

111 

$S1 

USER 

PT 

RCD 

273 

0 

0 

0 

149.4800 

-.9436 

16.4880 

111 

111 

$S1 

USER 

PT 

RCD 

274 

0 

0 

0 

149.4800 

24.1245 

-4.9995 

111 

111 

$S1 

USER 

PT 

RCD 

275 

0 

0 

0 

149.4800 

22.0625 

-8.9998 

111 

111 

$S1 

USER 

PT 

RCD 

276 

0 

0 

0 

149.4800 

19.2328 

-11.9998 

111 

111 

$S1 

USER 

PT 

RCD 

277 

0 

0 

0 

149.4800 

15.0000 

-14.4343 

111 

111 

$S1 

USER 

PT 

RCD 

278 

0 

0 

0 

149.4800 

9.3972 

-15.5198 

111 

111 

$S1 

USER 

PT 

RCD 

279 

0 

0 

0 

149.4800 

3.8030 

-16.0437 

111 

111 

$S1 

USER 

PT 

RCD 

280 

0 

0 

0 

149.4800 

1.5600 

-16.3084 

111 

111 

$S1 

USER 

PT 

RCD 

281 

0 

0 

0 

149.4800 

-.9436 

-16.4880 

111 

111 

$S1 

USER 

PT 

RCD 

282 

0 

0 

0 

162.7800 

20.8424 

.0000 

111 

111 

$S1 

USER 

PT 

RCD 

283 

0 

0 

0 

162.7800 

20.1915 

4.2271 

111 

111 

$S1 

USER 

PT 

RCD 

284 

0 

0 

0 

162.7800 

18.6665 

7.5112 

111 

111 

$S1 

USER 

PT 

RCD 

285 

0 

0 

0 

162.7800 

16.3676 

10.2830 

111 

111 

$S1 

USER 

PT 

RCD 

286 

0 

0 

0 

162.7800 

12.8861 

12.6611 

111 

111 

$S1 

USER 

PT 

RCD 

287 

0 

0 

0 

162.7800 

8.1095 

13.9902 

111 

111 

$S1 

USER 

PT 

RCD 

288 

0 

0 

0 

162.7800 

1.7300 

14.5972 

111 

111 

$S1 

USER 

PT 

RCD 

289 

0 

0 

0 

162.7800 

- .7750 

14.7893 

111 

111 

$S1 

USER 

PT 

RCD 

290 

0 

0 

0 

162.7800 

20.1915 

-4.2271 

111 

111 

$S1 

USER 

PT 

RCD 

291 

0 

0 

0 

162.7800 

18.6665 

-7.5112 

111 

111 

$S1 

USER 

PT 

RCD 

292 

0 

0 

0 

162.7800 

16.3676 

-10.2830 

111 

111 

$S1 

USER 

PT 

RCD 

293 

0 

0 

0 

162.7800 

12.8861 

-12.6611 

111 

111 

$S1 

USER 

PT 

RCD 

294 

0 

0 

0 

162.7800 

8.1095 

-13.9902 

111 

111 

$S1 

USER 

PT 

RCD 

295 

0 

0 

0 

162.7800 

1.7300 

-14.5972 

111 

111 

$S1 

USER 

PT 

RCD 

296 

0 

0 

0 

162.7800 

- .7750 

-14.7893 

111 

111 

$S1 

USER 

PT 

RCD 

297 

0 

0 

0 

176.0900 

16.2070 

.0000 

111 

111 

$S1 

USER 

PT 

RCD 

298 

0 

0 

0 

176.0900 

15.7459 

3.3538 

111 

111 

$S1 

USER 

PT 

RCD 

299 

0 

0 

0 

176.0900 

14.7690 

5.8030 

111 

111 

$S1 

USER 

PT 

RCD 

300 

0 

0 

0 

176.0900 

13.0393 

8.2886 

111 

111 

$S1 

USER 

PT 

RCD 

301 

0 

0 

0 

176.0900 

10.3600 

10.5422 

111 

111 

$S1 

USER 

PT 

RCD 

302 

0 

0 

0 

176.0900 

6.4710 

12.1328 

111 

111 

$S1 

USER 

PT 

RCD 

303 

0 

0 

0 

176.0900 

1.9200 

12.6490 

111 

111 

$S1 

USER 

PT 

RCD 

304 

0 

0 

0 

176.0900 

-.5825 

12.8499 

111 

111 

$S1 

USER 

PT 

RCD 

305 

0 

0 

0 

176.0900 

15.7459 

-3.3538 

111 

111 

$S1 

USER 

PT 

RCD 

306 

0 

0 

0 

176.0900 

14.7690 

-5.8030 

111 

111 

$S1 

USER 

PT 

RCD 

307 

0 

0 

0 

176.0900 

13.0393 

-8.2886 

111 

111 

$S1 

USER 

PT 

RCD 

308 

0 

0 

0 

176.0900 

10.3600 

-10.5422 

111 

111 

$S1 

USER 

PT 

RCD 

309 

0 

0 

0 

176.0900 

6.4710 

-12.1328 

111 

111 

$S1 

USER 

PT 

RCD 

310 

0 

0 

0 

176.0900 

1.9200 

-12.6490 

111 

111 

$S1 

USER 

PT 

RCD 

311 

0 

0 

0 

176.0900 

- .5825 

-12.8499 

111 

111 

$S1 

USER 

PT 

RCD 

312 

0 

0 

0 

189.3900 

11.4580 

.0000 

111 

111 

$S1 

USER 

PT 

RCD 

313 

0 

0 

0 

189.3900 

11.1631 

2.4536 

111 

111 

$S1 

USER 

PT 

RCD 

314 

0 

0 

0 

189.3900 

10.6543 

3.9995 

111 

111 

$S1 

USER 

PT 

RCD 

315 

0 

0 

0 

189.3900 

9.5424 

5.9995 

111 

111 

$S1 

USER 

PT 

RCD 

316 

0 

0 

0 

189.3900 

7.4991 

8.1426 

111 

111 

$S1 

USER 

PT 

RCD 

317 

0 

0 

0 

189.3900 

4.4203 

9.8081 

111 

111 

$S1 

USER 

PT 

RCD 

318 

0 

0 

0 

189.3900 

2.1100 

10.3972 

111 

111 

$S1 

USER 

PT 

RCD 

319 

0 

0 

0 

189.3900 

- .3881 

10.5117 

111 

111 

$S1 

USER 

PT 

RCD 

143 


320 

0 

0 

0 

189.3900 

11.1631 

-2.4536 

111 

111 

$S1 

USER 

PT 

RCD 

321 

0 

0 

0 

189.3900 

10.6543 

-3.9995 

111 

111 

$S1 

USER 

PT 

RCD 

322 

0 

0 

0 

189.3900 

9 . 5424 

-5.9995 

111 

111 

$S1 

USER 

PT 

RCD 

323 

0 

0 

0 

189.3900 

7.4991 

-8.1426 

111 

111 

$S1 

USER 

PT 

RCD 

324 

0 

0 

0 

189.3900 

4.4203 

-9.8081 

111 

111 

$S1 

USER 

PT 

RCD 

325 

0 

0 

0 

189.3900 

2.1100 

-10.3972 

111 

111 

$S1 

USER 

PT 

RCD 

326 

0 

0 

0 

189.3900 

-.3881 

-10.5117 

111 

111 

$S1 

USER 

PT 

RCD 

327 

0 

0 

0 

196.0500 

9.0705 

.0000 

111 

111 

$S1 

USER 

PT 

RCD 

328 

0 

0 

0 

196.0500 

8.8519 

2.0000 

111 

111 

$S1 

USER 

PT 

RCD 

329 

0 

0 

0 

196.0500 

8.0001 

4.3208 

111 

111 

$S1 

USER 

PT 

RCD 

330 

0 

0 

0 

196.0500 

6.0001 

6.8855 

111 

111 

$S1 

USER 

PT 

RCD 

331 

0 

0 

0 

196.0500 

3.7709 

8.3677 

111 

111 

$S1 

USER 

PT 

RCD 

332 

0 

0 

0 

196.0500 

2.1900 

9 . 0000 

111 

111 

$S1 

USER 

PT 

RCD 

333 

0 

0 

0 

196.0500 

-.3097 

9.2942 

111 

111 

$S1 

USER 

PT 

RCD 

334 

0 

0 

0 

196.0500 

8.8519 

-2.0000 

111 

111 

$S1 

USER 

PT 

RCD 

335 

0 

0 

0 

196.0500 

8.0001 

-4.3208 

111 

111 

$S1 

USER 

PT 

RCD 

336 

0 

0 

0 

196.0500 

6.0001 

-6.8855 

111 

111 

$S1 

USER 

PT 

RCD 

337 

0 

0 

0 

196.0500 

3.7709 

-8.3677 

111 

111 

$S1 

USER 

PT 

RCD 

338 

0 

0 

0 

196.0500 

2.1900 

-9.0000 

111 

111 

$S1 

USER 

PT 

RCD 

339 

0 

0 

0 

196.0500 

-.3097 

-9.2942 

111 

111 

$S1 

USER 

PT 

RCD 

340 

0 

0 

0 

202.7000 

6.7899 

.0000 

111 

111 

$S1 

USER 

PT 

RCD 

341 

0 

0 

0 

202.7000 

6.1991 

2.9995 

111 

111 

$S1 

USER 

PT 

RCD 

342 

0 

0 

0 

202 . 7000 

5.0311 

4.9995 

111 

111 

$S1 

USER 

PT 

RCD 

343 

0 

0 

0 

202.7000 

2.2800 

7 . 3000 

111 

111 

$S1 

USER 

PT 

RCD 

344 

0 

0 

0 

202.7000 

-.2170 

7.8560 

111 

111 

$S1 

USER 

PT 

RCD 

345 

0 

0 

0 

202.7000 

6.1991 

-2.9995 

111 

111 

$S1 

USER 

PT 

RCD 

346 

0 

0 

0 

202.7000 

5.0311 

-4.9995 

111 

111 

$S1 

USER 

PT 

RCD 

347 

0 

0 

0 

202.7000 

2.2800 

-7.3000 

111 

111 

$S1 

USER 

PT 

RCD 

348 

0 

0 

0 

202.7000 

-.2170 

-7.8560 

111 

111 

$S1 

USER 

PT 

RCD 

349 

0 

0 

0 

209.3600 

4.3978 

.0000 

111 

111 

$S1 

USER 

PT 

RCD 

350 

0 

0 

0 

209.3600 

4.0000 

2.4407 

111 

111 

$S1 

USER 

PT 

RCD 

351 

0 

0 

0 

209.3600 

2.3261 

5.0000 

111 

111 

$S1 

USER 

PT 

RCD 

352 

0 

0 

0 

209.3600 

-.1730 

6.2522 

111 

111 

$S1 

USER 

PT 

RCD 

353 

0 

0 

0 

209.3600 

4 . 0000 

-2.4407 

111 

111 

$S1 

USER 

PT 

RCD 

354 

0 

0 

0 

209.3600 

2.3261 

-5.0000 

111 

111 

$S1 

USER 

PT 

RCD 

355 

0 

0 

0 

209.3600 

-.1730 

-6.2522 

111 

111 

$S1 

USER 

PT 

RCD 

356 

0 

0 

0 

212.7610 

3.0416 

.0000 

111 

111 

$S1 

USER 

PT 

RCD 

357 

0 

0 

0 

215.5260 

1.3518 

.0000 

111 

111 

$S1 

USER 

PT 

RCD 

358 

0 

0 

0 

215.5260 

1.3518 

2.9995 

111 

111 

$S1 

USER 

PT 

RCD 

359 

0 

0 

0 

215.5260 

- .1112 

3.9998 

111 

111 

$S1 

USER 

PT 

RCD 

360 

0 

0 

0 

215.5260 

1.3518 

-2.9995 

111 

111 

$S1 

USER 

PT 

RCD 

361 

0 

0 

0 

215.5260 

-.1112 

-3.9998 

111 

111 

$S1 

USER 

PT 

RCD 

362 

0 

0 

0 

216.0100 

.0000 

.0000 

111 

111 

$S1 

USER 

PT 

RCD 

363 

0 

0 

0 

0.0000 

-0.5000 

0.0000 

001 

000 

$S1 

USER 

PT 

RCD 

364 

0 

0 

0 

1.7067 

-0.5000 

3.4593 

001 

000 

$S1 

USER 

PT 

RCD 

365 

0 

0 

0 

1.7067 

-0.5000 

-3.4593 

001 

000 

$S1 

USER 

PT 

RCD 

366 

0 

0 

0 

5.8027 

-0.6554 

6.9066 

001 

000 

$S1 

USER 

PT 

RCD 

367 

0 

0 

0 

5.8027 

-0.6554 

-6.9066 

001 

000 

$S1 

USER 

PT 

RCD 

368 

0 

0 

0 

11.6267 

-0.8141 

9.8173 

001 

000 

$S1 

USER 

PT 

RCD 

369 

0 

0 

0 

11.6267 

-0.8141 

-9.8173 

001 

000 

$S1 

USER 

PT 

RCD 

370 

0 

0 

0 

17.4507 

-0.9632 

11.9202 

001 

000 

$S1 

USER 

PT 

RCD 

371 

0 

0 

0 

17.4507 

-0.9632 

-11.9202 

001 

000 

$S1 

USER 

PT 

RCD 

372 

0 

0 

0 

23.2533 

-1.0935 

13.4892 

001 

000 

$S1 

USER 

PT 

RCD 

373 

0 

0 

0 

23.2533 

-1.0935 

-13.4892 

001 

000 

$S1 

USER 

PT 

RCD 

144 


374 

0 

0 

0 

34.9012 

-1.3093 

15.7464 

001 

000 

$S1 

USER 

PT 

ROD 

375 

0 

0 

0 

34.9012 

-1.3093 

-15.7464 

001 

000 

$S1 

USER 

PT 

ROD 

376 

0 

0 

0 

46.5493 

-1.4704 

17.2388 

001 

000 

$S1 

USER 

PT 

ROD 

377 

0 

0 

0 

46.5493 

-1.4704 

-17.2388 

001 

000 

$S1 

USER 

PT 

ROD 

378 

0 

0 

0 

58.1759 

-1.5806 

18.1876 

001 

000 

$S1 

USER 

PT 

RCD 

379 

0 

0 

0 

58.1759 

-1.5806 

-18.1876 

001 

000 

$S1 

USER 

PT 

RCD 

380 

0 

0 

0 

69.8239 

-1.6476 

18.7409 

001 

000 

$S1 

USER 

PT 

RCD 

381 

0 

0 

0 

69.8239 

-1.6476 

-18.7409 

001 

000 

$S1 

USER 

PT 

RCD 

382 

0 

0 

0 

75.6373 

-1.6636 

18.8704 

001 

000 

$S1 

USER 

PT 

RCD 

383 

0 

0 

0 

75.6373 

-1.6636 

-18.8704 

001 

000 

$S1 

USER 

PT 

RCD 

384 

0 

0 

0 

81.4506 

-1.6797 

18.9999 

001 

000 

$S1 

USER 

PT 

RCD 

385 

0 

0 

0 

81.4506 

-1.6797 

-18.9999 

001 

000 

$S1 

USER 

PT 

RCD 

386 

0 

0 

0 

90.8613 

-1.6880 

19.0881 

001 

000 

$S1 

USER 

PT 

RCD 

387 

0 

0 

0 

90.8613 

-1.6880 

-19.0881 

001 

000 

$S1 

USER 

PT 

RCD 

388 

0 

0 

0 

98.5000 

-1.6880 

19.0881 

001 

000 

$S1 

USER 

PT 

RCD 

389 

0 

0 

0 

98.5000 

-1.6880 

-19.0881 

001 

000 

$S1 

USER 

PT 

RCD 

390 

0 

0 

0 

103.8900 

-1.6880 

19.0881 

001 

000 

$S1 

USER 

PT 

RCD 

391 

0 

0 

0 

103.8900 

-1.6880 

-19.0881 

001 

000 

$S1 

USER 

PT 

RCD 

392 

0 

0 

0 

113.3800 

-1.6686 

18.8819 

001 

000 

$S1 

USER 

PT 

RCD 

393 

0 

0 

0 

113.3800 

-1.6686 

-18.8819 

001 

000 

$S1 

USER 

PT 

RCD 

394 

0 

0 

0 

122.8700 

-1.6492 

18.6758 

001 

000 

$S1 

USER 

PT 

RCD 

395 

0 

0 

0 

122.8700 

-1.6492 

-18.6758 

001 

000 

$S1 

USER 

PT 

RCD 

396 

0 

0 

0 

136.1740 

-1.5464 

17.8586 

001 

000 

$S1 

USER 

PT 

RCD 

397 

0 

0 

0 

136.1740 

-1.5464 

-17.8586 

001 

000 

$S1 

USER 

PT 

RCD 

398 

0 

0 

0 

149.4800 

-1.4436 

16.4880 

001 

000 

$S1 

USER 

PT 

RCD 

399 

0 

0 

0 

149.4800 

-1.4436 

-16.4880 

001 

000 

$S1 

USER 

PT 

RCD 

400 

0 

0 

0 

162.7800 

-1.2750 

14.7893 

001 

000 

$S1 

USER 

PT 

RCD 

401 

0 

0 

0 

162.7800 

-1.2750 

-14.7893 

001 

000 

$S1 

USER 

PT 

RCD 

402 

0 

0 

0 

176.0900 

-1.0825 

12.8499 

001 

000 

$S1 

USER 

PT 

RCD 

403 

0 

0 

0 

176.0900 

-1.0825 

-12.8499 

001 

000 

$S1 

USER 

PT 

RCD 

404 

0 

0 

0 

189.3900 

-0.8881 

10.5117 

001 

000 

$S1 

USER 

PT 

RCD 

405 

0 

0 

0 

189.3900 

-0.8881 

-10.5117 

001 

000 

$S1 

USER 

PT 

RCD 

406 

0 

0 

0 

196.0500 

-0.8097 

9.2942 

001 

000 

$S1 

USER 

PT 

RCD 

407 

0 

0 

0 

196.0500 

-0.8097 

-9.2942 

001 

000 

$S1 

USER 

PT 

RCD 

408 

0 

0 

0 

202.7000 

-0.7170 

7.8560 

001 

000 

$S1 

USER 

PT 

RCD 

409 

0 

0 

0 

202.7000 

-0.7170 

-7.8560 

001 

000 

$S1 

USER 

PT 

RCD 

410 

0| 

0 

0 

209.3600 

-0.6730 

6.2522 

001 

000 

$S1 

USER 

PT 

RCD 

411 

0 

0 

0 

209.3600 

-0.6730 

-6.2522 

001 

000 

$S1 

USER 

PT 

RCD 

412 

0 

0 

0 

215.5260 

-0.6112 

3.9998 

001 

000 

$S1 

USER 

PT 

RCD 

413 

0 

0 

0 

215.5260 

-0.6112 

-3.9998 

001 

000 

$S1 

USER 

PT 

RCD 

414 

0 

0 

0 

216.0100 

-0.5000 

0.0000 

001 

000 

$S1 

USER 

PT 

RCD 

415 

0 

0 

0 

0.0000 

-0.8770 

0 . 0000 

000 

000 

$S1 

USER 

PT 

RCD 

416 

0 

0 

0 

1.7067 

-0.8770 

3.4593 

000 

000 

$S1 

USER 

PT 

RCD 

417 

0 

0 

0 

1.7067 

-0.8770 

-3.4593 

000 

000 

$S1 

USER 

PT 

RCD 

418 

0 

0 

0 

5.8027 

-1.0324 

6.9066 

000 

000 

$S1 

USER 

PT 

RCD 

419 

0 

0 

0 

5.8027 

-1.0324 

-6.9066 

000 

000 

$S1 

USER 

PT 

RCD 

420 

0 

0 

0 

11.6267 

-1.1911 

9.8173 

000 

000 

$S1 

USER 

PT 

RCD 

421 

0 

0 

0 

11.6267 

-1.1911 

-9.8173 

000 

000 

$S1 

USER 

PT 

RCD 

422 

0 

0 

0 

17.4507 

-1.3402 

11.9202 

000 

000 

$S1 

USER 

PT 

RCD 

423 

0 

0 

0 

17.4507 

-1.3402 

-11.9202 

000 

000 

$S1 

USER 

PT 

RCD 

424 

0 

0 

0 

23.2533 

-1.4705 

13.4892 

000 

000 

$S1 

USER 

PT 

RCD 

425 

0 

0 

0 

23.2533 

-1.4705 

-13.4892 

000 

000 

$S1 

USER 

PT 

RCD 

426 

0 

0 

0 

34.9012 

-1.6863 

15.7464 

000 

000 

$S1 

USER 

PT 

RCD 

427 

0 

0 

0 

34.9012 

-1.6863 

-15.7464 

000 

000 

$S1 

USER 

PT 

RCD 

145 


428 

0 

0 

0 

46.5493 

-1.8474 

17.2388 

000 

000 

$S1 

USER 

FT 

ROD 

429 

0 

0 

0 

46.5493 

-1.8474  -17.2388 

000 

000 

$S1 

USER 

FT 

ROD 

430 

0 

0 

0 

58.1759 

-1.9576 

18.1876 

000 

000 

$S1 

USER 

FT 

ROD 

431 

0 

0 

0 

58.1759 

-1.9576  -18.1876 

000 

000 

$S1 

USER 

FT 

RCD 

432 

0 

0 

0 

69.8239 

-2.0246 

18 . 7409 

000 

000 

$S1 

USER 

FT 

ROD 

433 

0 

0 

0 

69.8239 

-2.0246  -18.7409 

000 

000 

$S1 

USER 

FT 

RCD 

434 

0 

0 

0 

75.6373 

-2.0406 

18.8704 

000 

000 

$S1 

USER 

FT 

RCD 

435 

0 

0 

0 

75.6373 

-2.0406  -18.8704 

000 

000 

$S1 

USER 

FT 

RCD 

436 

0 

0 

0 

81.4506 

-2.0567 

18.9999 

000 

000 

$S1 

USER 

FT 

RCD 

437 

0 

0 

0 

81.4506 

-2.0567  -18.9999 

000 

000 

$S1 

USER 

FT 

RCD 

438 

0 

0 

0 

90.8613 

-2.0650 

19.0881 

000 

000 

$S1 

USER 

FT 

RCD 

439 

0 

0 

0 

90.8613 

-2.0650  -19.0881 

000 

000 

$S1 

USER 

FT 

RCD 

440 

0 

0 

0 

98.5000 

-2.0650 

19.0881 

000 

000 

$S1 

USER 

FT 

RCD 

441 

0 

0 

0 

98.5000 

-2.0650  -19.0881 

000 

000 

$S1 

USER 

FT 

RCD 

442 

0 

0 

0 

103.8900 

-2.0650 

19.0881 

000 

000 

$S1 

USER 

FT 

RCD 

443 

0 

0 

0 

103.8900 

-2.0650  -19.0881 

000 

000 

$S1 

USER 

FT 

RCD 

444 

0 

0 

0 

113.3800 

-2.0456 

18.8819 

000 

000 

$S1 

USER 

FT 

RCD 

445 

0 

0 

0 

113.3800 

-2.0456  -18.8819 

000 

000 

$S1 

USER 

FT 

RCD 

446 

0 

0 

0 

122.8700 

-2.0262 

18.6758 

000 

000 

$S1 

USER 

FT 

RCD 

447 

0 

0 

0 

122.8700 

-2.0262  -18.6758 

000 

000 

$S1 

USER 

FT 

RCD 

448 

0 

0 

0 

136.1740 

-1.9234 

17.8586 

000 

000 

$S1 

USER 

FT 

RCD 

449 

0 

0 

0 

136.1740 

-1.9234  -17.8586 

000 

000 

$S1 

USER 

FT 

RCD 

450 

0 

0 

0 

149.4800 

-1.8206 

16.4880 

000 

000 

$S1 

USER 

FT 

RCD 

451 

0 

0 

0 

149.4800 

-1.8206  -16.4880 

000 

000 

$S1 

USER 

FT 

RCD 

452 

0 

0 

0 

162.7800 

-1.6520 

14.7893 

000 

000 

$S1 

USER 

FT 

RCD 

453 

0 

0 

0 

162.7800 

-1.6520  -14.7893 

000 

000 

$S1 

USER 

FT 

RCD 

454 

0 

0 

0 

176.0900 

-1.4595 

12.8499 

000 

000 

$S1 

USER 

FT 

RCD 

455 

0 

0 

0 

176.0900 

-1.4595  -12.8499 

000 

000 

$S1 

USER 

FT 

RCD 

456 

0 

0 

0 

189.3900 

-1.2651 

10.5117 

000 

000 

$S1 

USER 

FT 

RCD 

457 

0 

0 

0 

189.3900 

-1.2651  -10.5117 

000 

000 

$S1 

USER 

FT 

RCD 

458 

0 

0 

0 

196.0500 

-1.1867 

9.2942 

000 

000 

$S1 

USER 

FT 

RCD 

459 

0 

0 

0 

196.0500 

-1.1867 

-9.2942 

000 

000 

$S1 

USER 

FT 

RCD 

460 

0 

0 

0 

202.7000 

-1.0940 

7.8560 

000 

000 

$S1 

USER 

FT 

RCD 

461 

0 

0 

0 

202 . 7000 

-1.0940 

-7.8560 

000 

000 

$S1 

USER 

FT 

RCD 

462 

0 

0 

0 

209.3600 

-1.0500 

6.2522 

000 

000 

$S1 

USER 

FT 

RCD 

463 

0 

0 

0 

209.3600 

-1.0500 

-6.2522 

000 

000 

$S1 

USER 

FT 

RCD 

464 

0 

0 

0 

215.5260 

-0.9882 

3.9998 

000 

000 

$S1 

USER 

FT 

RCD 

465 

0 

0 

0 

215.5260 

-0.9882 

-3.9998 

000 

000 

$S1 

USER 

FT 

RCD 

466 

0 

0 

0 

216.0100 

-0.8770 

0 . 0000 

000 

000 

$S1 

USER 

FT 

RCD 

363 

415 

365 

210 

1 

0.0 

0.0 

0.0 

1 

0 

$T2 

BEAM 

ELEM 

364 

416 

363 

210 

1 

0.0 

0.0 

0.0 

1 

0 

$T2 

BEAM 

ELEM 

366 

418 

364 

210 

1 

0.0 

0.0 

0.0 

1 

0 

$T2 

BEAM 

ELEM 

368 

420 

366 

210 

1 

0.0 

0.0 

0.0 

1 

0 

$T2 

BEAM 

ELEM 

370 

422 

368 

210 

1 

0.0 

0.0 

0.0 

1 

0 

$T2 

BEAM 

ELEM 

372 

424 

370 

210 

1 

0.0 

0.0 

0.0 

1 

0 

$T2 

BEAM 

ELEM 

374 

426 

372 

210 

1 

0.0 

0.0 

0.0 

1 

0 

$T2 

BEAM 

ELEM 

376 

428 

374 

210 

1 

0.0 

0.0 

0.0 

1 

0 

$T2 

BEAM 

ELEM 

378 

430 

376 

210 

1 

0.0 

0.0 

0.0 

1 

0 

$T2 

BEAM 

ELEM 

380 

432 

378 

210 

1 

0.0 

0.0 

0.0 

1 

0 

$T2 

BEAM 

ELEM 

382 

434 

380 

210 

1 

0.0 

0.0 

0.0 

1 

0 

$T2 

BEAM 

ELEM 

384 

436 

382 

210 

1 

0.0 

0.0 

0.0 

1 

0 

$T2 

BEAM 

ELEM 

386 

438 

384 

210 

1 

0.0 

0.0 

0.0 

1 

0 

$T2 

BEAM 

ELEM 

388 

440 

386 

210 

1 

0.0 

0.0 

0.0 

1 

0 

$T2 

BEAM 

ELEM 

390 

442 

388 

210 

1 

0.0 

0.0 

0.0 

1 

0 

$T2 

BEAM 

ELEM 

146 


392  444  390  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

394  446  392  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

396  448  394  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

398  450  396  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

400  452  398  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

402  454  400  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

404  456  402  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

406  458  404  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

408  460  406  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

410  462  408  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

412  464  410  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

414  466  412  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

413  465  414  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

411  463  413  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

409  461  411  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

407  459  409  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

405  457  407  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

403  455  405  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

401  453  403  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

399  451  401  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

397  449  399  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

395  447  397  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

393  445  395  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

391  443  393  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

389  441  391  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

387  439  389  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

385  437  387  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

383  435  385  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

381  433  383  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

379  431  381  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

377  429  379  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

375  427  377  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

373  425  375  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

371  423  373  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

369  421  371  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

367  419  369  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

365  417  367  210  1  0.0  0.0  0.0  1  0  $T2  BEAM  ELEM 

1  3  2  320  26  $T3  TRI  ELEM 

1  4  3  320  26  $T3  TRI  ELEM 

1  2  5  320  26  $T3  TRI  ELEM 

1  5  6  320  26  $T3  TRI  ELEM 

2  8  7  320  15  $T3  TRI  ELEM 

2  9  8  320  15  $T3  TRI  ELEM 

2  3  9  320  15  $T3  TRI  ELEM 

3  10  9  320  15  $T3  TRI  ELEM 

2  7  12  320  15  $T3  TRI  ELEM 

2  12  13  320  15  $T3  TRI  ELEM 

2  13  5  320  15  $T3  TRI  ELEM 

5  13  14  320  15  $T3  TRI  ELEM 

7  17  16  320  15  $T3  TRI  ELEM 

9  20  19  320  15  $T3  TRI  ELEM 

10  20  9  320  15  $T3  TRI  ELEM 

10  21  20  320  15  $T3  TRI  ELEM 

7  16  23  320  15  $T3  TRI  ELEM 
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Appendix  H:  UPRESS  Subroutine  for  Finite  Element 
Models  6  and  7 


The  following  subroutine  was  used  for  applying  the  pressure 
loading  in  the  Model  6  and  7  finite  element  runs .  The  subroutine  is 
written  using  guidelines  in  the  STAGS  users  manual  [16] .  The 
subroutine  is  then  compiled  and  linked  to  the  main  program  before 
execution. 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


20 

30 


C 


SUBROUTINE  UPRESS <T , PA , PB , lUNIT , lELT , X , Y , Z , LIVE , PRESS ) 


Pressure  distribution  subroutine  for  the  modified 
boundary  conditions.  Uses  element  number  to  determine 
element  pressure  loading.  Modified  from  Models  1-5 
because  of  the  addition  of  the  52  rubber  beam  elements 
(not  loaded)  in  models  6  and  7.  The  52  beams  are 
elements  1  through  52.  The  next  362  elements  correspond 
to  the  original  finite  element  model.  All  elements 
after  this  are  plate  elements  used  to  model  the  aluminum 
plate  in  the  modified  boundary  conditions. 

Needs  input  file  'loads.dat'  which  contains  the  twenty 
loads  written  one  on  each  line. 


INTEGER  I.NIELT 
REAL  C(20) ,PA 

OPEN  (UNIT-l, FILE-' LOADS. DAT' , STATUS- ' OLD ' ) 
LIVE-1 

DO  20  1-1,20 
READ(1,30)  C(I) 

CONTINUE 
FORMAT (FI 5. 7) 

REWIND  (UNIT-1) 

IF  ((IELT.LE.52) .OR. (IELT.GE.415))  THEN 
PRESS-0.0 


ELSE 

NIELT-IELT-52 
IF  ((NIELT.LE. 
+.0R. (NIELT.GE. 
+.0R. (NIELT.EQ. 
+.0R. (NIELT.EQ. 
+.0R. (NIELT.GE. 
+.OR. (NIELT.GE. 
+.0R. (NIELT.EQ. 
+.0R. (NIELT.GE. 


2) 

5. AND. NIELT.LE.  8) 
55) 

13) 

57. AND. NIELT.LE.  59) 
14. AND. NIELT.LE.  16) 
21) 

63. AND. NIELT.LE.  67) 


156 


+.0R. (NIELT.GE.  22 . AND.NIELT.LE.  23) 

+.0R. (NIELT.GE.  73. AND.NIELT.LE.  78) 

+. OR. (NIELT.GE.  27 .AND.NIELT.LE.  29) 

+. OR. (NIELT.GE.  85 .AND.NIELT.LE.  92))  PRESS-C(1)*PA 
C 

IF  ( (NIELT . GE . 101 .AND.NIELT.LE. 104) 

+.0R. (NIELT.GE. 117. AND.NIELT.LE. 120) 

+ . OR . (NIELT . GE . 133 . AND . NIELT . LE . 136 ) )  PRESS-C ( 2 ) *PA 
C 

IF  ((NIELT.GE. 105. AND.NIELT.LE. 108) 

+.OR. (NIELT.GE. 121. AND.NIELT.LE. 124) 

+.OR. (NIELT.GE. 137. AND.NIELT.LE. 140))  PRESS-C ( 3 )*PA 
C 

IF  ( (NIELT . GE . 149 . AND . NIELT . LE . 152 ) 

+.OR. (NIELT.GE. 157. AND.NIELT.LE. 160) 

+.OR. (NIELT.GE. 181. AND.NIELT.LE. 184) 

+.OR. (NIELT.GE. 197. AND.NIELT.LE. 200) 

+.OR. (NIELT.GE. 205. AND.NIELT.LE. 208))  PRESS-C(4)*PA 
C 

IF  ((NIELT.GE. 153. AND.NIELT.LE. 156) 

+.OR. (NIELT.GE. 161. AND.NIELT.LE. 164) 

+.OR. (NIELT.GE. 185. AND.NIELT.LE. 188) 

+.OR. (NIELT.GE. 201. AND.NIELT.LE. 204) 

+.OR. (NIELT.GE. 209. AND.NIELT.LE. 212))  PRESS-C ( 5 )*PA 
C 

IF  ((NIELT.GE. 229. AND.NIELT.LE. 232) 

+.OR. (NIELT.GE. 237. AND.NIELT.LE. 240) 

+. OR . (NIELT . GE . 261 . AND . NIELT . LE . 264) )  PRESS-C ( 6 ) *PA 
C 

IF  ((NIELT.GE. 233. AND.NIELT.LE. 236) 

+.OR. (NIELT.GE. 241. AND.NIELT.LE. 244) 

+.OR. (NIELT.GE. 265. AND.NIELT.LE. 268))  PRESS-C ( 7) *PA 
C 

IF  ((NIELT.GE. 277. AND.NIELT.LE. 281) 

+.OR. (NIELT.GE. 293. AND.NIELT.LE. 297) 

+ . OR . (NIELT . GE . 307 . AND . NIELT . LE . 310) )  PRESS-C ( 8 ) *PA 
C 

IF  ((NIELT.GE. 282. AND.NIELT.LE. 284) 

+.OR. (NIELT.GE. 298. AND.NIELT.LE. 299) 

+. OR. (NIELT. EQ.  33) 

+.OR. (NIELT.GE. 311. AND.NIELT.LE. 313))  PRESS-C (9 )*PA 
C 

IF  ((NIELT.GE. 321. AND.NIELT.LE. 327) 

+.OR. (NIELT.GE. 335. AND.NIELT.LE. 337) 

+. OR. (NIELT. EQ.  35) 

+.OR. (NIELT.GE. 347. AND.NIELT.LE. 350) 

+. OR. (NIELT.GE.  37 .AND.NIELT.LE.  38) 

+.OR. (NIELT.GE. 355. AND.NIELT.LE. 357) 

+. OR. (NIELT. EQ.  41) 

+. OR. (NIELT.GE.  43 .AND . NIELT . LE .  45) 

+.OR. (NIELT. EQ. 361) 

+. OR. (NIELT. EQ.  49) 

+. OR. (NIELT.GE.  51 .AND.NIELT.LE.  52))  PRESS-C(10)*PA 
C 
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IF  ((NIELT.GE.  3 .AND.NIELT.LE.  4) 

+.0R. (NIELT.GE.  9 .AND.NIELT.LE.  12) 

+.0R. (NIELT.EQ.  56) 

+.0R. (NIELT.EQ.  17) 

+. OR. (NIELT.GE.  60. AND.NIELT.LE.  62) 

+. OR. (NIELT.GE.  18. AND.NIELT.LE.  20) 

+. OR. (NIELT.EQ.  24) 

+. OR. (NIELT.GE.  68. AND.NIELT.LE.  72) 

+. OR. (NIELT.GE.  25 .AND.NIELT.LE.  26) 

+. OR. (NIELT.GE.  79 .AND.NIELT.LE.  84) 

+. OR. (NIELT.GE.  30. AND.NIELT.LE.  32) 

+. OR. (NIELT.GE.  93 .AND.NIELT.LE. 100) )  PRESS-C(11)*PA 
C 

IF  ((NIELT.GE. 109. AND.NIELT.LE. 112) 

+.OR. (NIELT.GE. 125. AND.NIELT.LE. 128) 

+ . OR . (NIELT . GE . 141 . AND . NIELT . LE . 144) )  PRESS-C ( 12 ) *PA 
C 

IF  ((NIELT.GE. 113. AND.NIELT.LE. 116) 

+.OR. (NIELT.GE. 129. AND.NIELT.LE. 132) 

+  . OR .( NIELT . GE . 145 . AND . NIELT . LE . 148 ) )  PRES S-C ( 1 3 ) *PA 
C 

IF  ((NIELT.GE. 165. AND.NIELT.LE. 168) 

+.OR. (NIELT.GE. 173. AND.NIELT.LE. 176) 

+.OR. (NIELT.GE. 189. AND.NIELT.LE. 192) 

+.OR. (NIELT.GE. 213. AND.NIELT.LE. 216) 

+.OR. (NIELT.GE. 221. AND.NIELT.LE. 224))  PRESS-C(14)*PA 
C 

IF  ((NIELT.GE. 169. AND.NIELT.LE. 172) 

+.OR. (NIELT.GE. 177. AND.NIELT.LE. 180) 

+.OR. (NIELT.GE. 193. AND.NIELT.LE. 196) 

+.OR. (NIELT.GE. 217. AND.NIELT.LE. 220) 

+.OR. (NIELT.GE. 225. AND.NIELT.LE. 228))  PRESS-C ( 15 )*PA 
C 

IF  ( (NIELT . GE . 245 . AND . NIELT . LE . 248 ) 

+,OR. (NIELT.GE. 253. AND.NIELT.LE. 256) 

+.OR. (NIELT.GE. 269. AND.NIELT.LE. 272))  PRESS-C ( 16 )*PA 
C 

IF  ((NIELT.GE. 249. AND.NIELT.LE. 252) 

+.OR. (NIELT.GE. 257. AND.NIELT.LE. 260) 

+.OR. (NIELT.GE. 274. AND.NIELT.LE. 276))  PRESS-C ( 17 )*PA 
C 

IF  ((NIELT.GE. 285. AND.NIELT.LE. 289) 

+ . OR. (NIELT . GE . 300 . AND . NIELT . LE . 304) 

+ . OR . (NIELT . GE . 314 . AND . NIELT . LE . 317 ) )  PRESS=C ( 18 ) *PA 
C 

IF  ((NIELT.GE. 290. AND.NIELT.LE. 292) 

+.OR. (NIELT.GE. 305. AND.NIELT.LE. 306) 

+. OR. (NIELT.EQ.  34) 

+. OR . (NIELT . GE . 318 . AND . NIELT . LE . 320 ) )  PRESS-C (19) *PA 
C 

IF  ( (NIELT . GE . 328 .AND . NIELT . LE . 334) 

+.OR. (NIELT.GE. 341. AND.NIELT.LE. 346) 

+. OR. (NIELT.EQ.  36) 

+.OR. (NIELT.GE. 351. AND.NIELT.LE. 354) 
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+.0R. (NIELT.GE.  39 .AND.NIELT.LE.  40) 

+.OR. (NIELT.GE. 358. AND.NIELT.LE. 360) 

+.OR. (NIELT.EQ.  42) 

+. OR. (NIELT.GE.  46. AND.NIELT.LE.  48) 

+.OR. (NIELT.EQ. 362) 

+. OR. (NIELT.EQ.  50) 

+. OR. (NIELT.GE.  53 .AND.NIELT.LE.  54))  PRESS-C(20)*PA 

ENDIF 

RETURN 

END 
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